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Abstract: Spectra of high-symmetry molecules contain fine and superfine level cluster 
structure related to J-tunneling between hills and valleys on rovibronic energy surfaces 
(RES). Such graphic visualizations help disentangle multi-level dynamics, selection rules, 
and state mixing effects including widespread violation of nuclear spin symmetry species. 
A review of RES analysis compares it to that of potential energy surfaces (PES) used 
in Born-Oppenheimer approximations. Both take advantage of adiabatic coupling in 
order to visualize Hamiltonian eigensolutions. RES of symmetric and D2 asymmetric top 
rank-2-tensor Hamiltonians are compared with Oh spherical top rank-4-tensor fine-structure 
clusters of 6-fold and 8-fold tunneling multiplets. Then extreme 12-fold and 24-fold 
multiplets are analyzed by RES plots of higher rank tensor Hamiltonians. Such extreme 
clustering is rare in fundamental bands but prevalent in hot bands, and analysis of its 
superfine structure requires more efficient labeling and a more powerful group theory. This 
is introduced using elementary examples involving two groups of order-6 {Cq and i^s^Cs^), 
then applied to families of Oh clusters in SFq spectra and to extreme clusters. 

Keywords: symmetry; molecular dynamics; tunnehng; level clusters 



1. Overview of Eigensolution Techniques for Symmetric Molecules 

A key mathematical technique for atomic or molecular physics and quantum chemistry is matrix 
diagonalization for quantum eigensolution. As computers become faster and more available, more 
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problems of chemical physics are framed in terms of choosing bases for eigensolution of time evolution 
operators or Hamiltonian generator matrices. The resulting eigenvectors and eigenvalues are Fourier 
amplitudes and frequencies that combine to give all possible dynamics in a given basis choice. 

Despite the increasing utility and power of computer diagonalization, it remains a "black box" of 
processes quite unlike the complex natural selection by wave interference that we imagine nature uses to 
arrive at its quantum states. Diagonalization uses numerical tricks to reduce each A^-by-A^ matrix to 
values and stationary eigenstates, but the artificial processes may seem as opaque as nature itself with 
little or no physical insight provided by A^^ — A^ eigenvector components. We are thus motivated to seek 
ways to visualize more of the physics of molecular eigensolutions and their spectra. This leads one to 
explore digital graphical visualization techniques that provide insight as well as increased computational 
power and thereby complement numerically intensive approaches [1]. 

Before describing tensor eigensolution techniques and rovibronic energy surfaces (RES), a brief 
review is given of potential energy surface (PES) to put the tensor RES in a historical and methodological 
context. This includes some background on semiclassical approximations of tensor algebra that help 
explain rotational level clustering and are used to develop the RES graphical tools. Section 2 reviews 
how RES apply to symmetric and asymmetric top molecules. This serves to motivate the application 
of RES to more complicated molecules of higher symmetry. Section 3 contains a graphical analysis 
of octahedral RES and an introductory review of level clusters (fine structure) having 6-fold and 8-fold 
quasi-degeneracy (superfine structure) due to rank-4 tensor Hamiltonians. Following this is a discussion 
of mixed-rank tensors that exhibit 12-fold and 24-fold monster-clusters. The latter have only recently 
been seen in highly excited rovibrational spectra [2] and present challenging problems of symmetry 
analysis to sort out a plethora of tunneling resonances and parameters for so many resonant states. 

Following introductory Section 4, these problems are addressed in Sections 6-8 by redeveloping 
group algebraic symmetry analysis into a more physically direct and elegantly powerful approach. It 
uses underlying duality between internal and external symmetry states and their operations. Duality 
is introduced using the simplest order-6 symmetry groups Cq and D^^Csy before applying it to Oh 
symmetric monster-clusters in Section 8. Monsters in REES-polyad bands are shown in final Section 9. 

The direct approach to symmetry starts by viewing a group product table as a Hamiltonian matrix H 
representing an H operator that is a linear combination of group operators g^^^ with a set of ortho-complete 
tunneling coefficients gk labeling each tunneling path. A main idea is that symmetry operators "know" 
the eigensolutions of their algebra and thus of all Hamiltonian and evolution operators made of g/e's. 

7.7. Computer Graphical Techniques 

Several graphical techniques and procedures exist for gaining spectral insight. One of the oldest is 
the Born-Oppenheimer approximate (BOA) potential energy surface (PES) that is a well-established 
tool for disentangling vibrational-electronic (vibronic) dynamics. While BOA-PES predate the digital 
age by decades, their calculation and display is made practical by computer. More recent are studies 
of phase portraits and wavepacket propagation techniques to follow high-z/ vibrational dynamics and 
chemical pathways for dissociation or re-association [3,4]. This includes BOA-breakdown states in 
which a system evolves on multiply interfering PES paths. Dynamic Jahn-Teller-Renner effects involve 
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multi-BOA-PES states in molecules and solids. Examples in recent works [5,6] include coherent 
photo-synthesis [7]. 

Visualizing eigensolutions and spectra in crystalline solids is helped by bands of dispersion functions 
in reciprocal frequency- versus-wavevector space. Fermi-sea contours are used to analyze de Haas-van 
Alphen effects and more recently in understanding quantum Hall effects. Analogy between band theory 
of solids and molecular rovibronic clusters is made in Section 4 and 7. 

Visualization of molecular rotational, rovibrational, and rovibronic eigensolutions and spectra is 
the subject of this work and involves the rotational energy surface (RES). As described below, an 
RES is a multipole expansion plot of an effective Hamiltonian in rotational momentum space. Ultra 
sensitivity of vibronic states to rotation lets the RES expose subtle and unexpected physics. Multi-RES 
or rovibronic energy eigenvalue surfaces (REES) have conical intersections analogous to Jahn-Teller 
PES (See Section 9). 

The RES was introduced about thirty years ago [8] to analyze spectral fine structure of high resolution 
spectral bands in molecules of high symmetry including PH^ [9] XDH^ and XD2H molecules [10], 
tetrahedral (P4) [H], tetrafluorides (CF4 and SiF^) [12], hexafluorides {SFq, Mo{CO)q and 
UFq) [13-16], cubane (C^H^), and buckyball (Ceo) [15,16] and predicted major mixing of Herzberg 
rovibronic species. Recently RES have been extended to help understand the dynamics and spectra of 
fluxional rotors [17] or "floppy" molecules such as methyl-complexes [18] and vibrational overtones [2]. 

Each of the techniques and particularly the RES -based ones described below depend upon the 
key wave functional properties of stationary phase, adiabatic invariance, and the spacetime symmetry 
underlying quantum theory. Additional symmetry (point group, space group, exchange, gauge, etc.) of a 
molecular system introduces additional resonance. Symmetry tends to make graphical techniques even 
more useful since they help clarify resonant phenomenal dynamics and symmetry labeling [19]. 

1.1.1. Vibronic Born-Openheimer Approximate Potential Energy Surfaces (BOA-PES) 

A BOA-PES depends on an adiabatic invariance of each electronic wavefunction to nuclear vibration. 
It is often said that the electrons are so much faster than nuclei that the system "sticks" to a particular 
PES that electrons provide. Perhaps a better criterion would be that the Fourier spectrum associated with 
nuclear motion does not overlap that of an electronic transition to another energy level. Nuclei often 
provide stable configurations that quantize electronic energy into levels separated by gaps much wider 
than that of low lying vibrational "phonon" states. 

A BOA wavefunction is a peculiarly entangled outer product = 77?/; of a nuclear factor wavefunction 
r7^(e) (X . . . ) whose quantum labels u{e) depend on electronic quantum numbers e = nlm, etc. while 
the electronic factor wave . . . ) ^ function whose electron coordinates • • • depend 

adiabatically on nuclear vibrational coordinates {X . . .) of PES V^{X . . .) for electron bond state e. 

iTf f ^electron -\r nuclei \ 

... A • • • j = 

V^e(^(x... )•••)• Vy{e) • • •) BOA-Eutauglcd Product (la) 
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^^electron -^nuclei 



. . . ) • r]jy{X . . .) Unentangled Product 



(lb) 



The adiabatic convenience of a single product Equation (la) with a vibration eigenfunction 
7y^(e)(X . . . ) on a single PES function Ve{X . . .) is welcome but comes at a price; a BOA-entangled 
coordinate- state is not a simple bra-ket wavefunction product Equation (lb) of bra-bra {x . . .\{X . . .\ 
position and ket-ket \ipe)\Viy) state. Symmetry operator product analysis of Equation (lb) is well known. 
Symmetry of Equation (la) depends on rotational B OA-relativity of its parts. Vibronic BOA-PES 
generalize to ro vibronic RES by accounting for rotational relations. 

1.1.2. Rovibronic BOA Rotational Energy Surfaces (BOA-RES) 

The rotational energy surface (RES) can be seen as a generalization of adiabatic BOA wave 
Equation (la) to Equation (2) below that includes rotational motion. Here one treats vibronic motion 
as having the "fast" degrees of freedom while rotational coordinates B (e.g., Euler angle {a/3j) for 
semi-rigid molecules) play the "slow" semi-classical role vis-a-vis the "faster" adiabatic vibration or 
vibronic states. 



In Equation (2), the wave factors of each motion are ordered fast-to- slow going left- to-right. As in 
Equation (la) each wave-factor quantum number depends on quanta in "faster" wave-factors written 
to its left, but each coordinate has adiabatic dependence on coordinates in "slower" factors written to 
its right. 

The Q in Equation (2) denotes vibrational normal coordinates (gi, ^2, • • • Qm) and u denotes their 
quanta (z^i, ^^2, • • • ^m)- The number m = 3A^ — 6 of modes of an A^-atom semi-rigid molecule has 
subtracted 3 translational and 3 rotational coordinates. Each mode qk assumes an adiabatic BOA 
dependency on overall translation and rotation B known as the Eckart conditions. (Here we will 
ignore translation.) 

RES are multipole expansion plots of effective BOA energy tensors for each quantum value of 
vibronic z/(e) and conserved total angular momentum J. Choices of effective energy tensors depend on 
the level of adiabatic approximation. So do the choices of spaces in which RES are plotted. Elementary 
examples of model BOA waves, tensors, and RES for rigid or semi-rigid molecules are discussed below. 

7.2. Lab-Frame Coupling vs. Body Frame Constriction 

Wave p^^Q^otation^^ j^j, ^ hdiVQ, rigid symmetric-top [ifj = 1 = r]) molecule is a Wigner L)^-function. 



UxiQ...s...) ■■■)■ •••[©•••])• PJKe)](©™* • • • ) 



(2) 



Pj(€>) = Pj,m,k("/57) = ^M,x("/57)Vnorm 



norm = [J] = 2J + 1 



(3) 



Total angular momentum J is J = i? for a bare rotor. Bare lab-frame z-component is labeled M = m. 
Its body-frame z-component is labeled K = M = n. m and n range from +R to — i? in integral steps. 
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Entangled BOA product Equation (2) mates vibronic factor Equation (la) with a rotor factor Pj = 
Pj,M,K Equation (3). Now J and K = M depend on total vibronic momentum / and its body z- 
component ft in '^^^{e) = 

^JHe)] = ^^L(e) • = % • Pj,M,K = % • ^mV^M (4) 

Disentangled product in Equation (lb) of lab-based vibronic wave and bare rotor p^^^ of 
Equation (3) is coupled by Clebsch-Gordan Coefficients C^^^ into a wave $}(^ of total J = R + I, 

+ / — 1,... or|i? — /| and M = /x + m by sum Equation (5) over lab z-angular bare rotor momenta 
m and lab vibronic /x bases. 

*M = E ■ = E ^i-M^i -^ZVW] (M = /I + m = const.) (5) 

A B OA-entangled wave in Equation (la) or Equation (4) requires more serious surgery in order to 
survive as a viable theoretical entity. BOA vibronic waves are not merely coupled as in Equation (5) 
to a rotor, they are adiabatically "glued" or constricted to the intrinsic molecular rotor frame. (A rotor 
is "BOA-constricted" by its vibronic wave much as a boa-constrictor rides its writhing prey as the two 
rotate together.) 

A remarkable property of quantum rotor operator algebra is that Wigner i^^-waves in Equation (3) are 
also transformation matrices that relate rotating body-fixed BOA ^!^j^{hody) into the lab-fixed "^^^{lah). 

¥-^{hody) = E K(^ab)^l^{a^^) (6a) 
Ki^cib) = E %{body)^l^a/3^) (6b) 

This rotational wave relativity is a subset of Lorentz-Einstein-Minkowski space-time-frame relativity 
that uses symmetry algebra to keep track of the invariant sub- spaces (eigensolutions). D-Matrices 
underlie all tensor operators, their eigenfunctions and their eigenvalues and are a non-Abelian 
(non-commutative) generalization of plane waves (i^*(r) = {r\k) = e^^^ underlying Fourier operator 
analysis. Details of this connection comprise the later Section 5. 

Of particular importance to RES theory is the Wigner-Eckart factorization lemma that relates 
Clebsch-Gordan C^^j^ to Wigner- i^'s and transforms coupled wave Equation (5) to BOA-constricted 
wave Equation (4). 

j d{a^^)&%{a^^)^^,{ah)^i,^{ah) = ^ffr^u^^i (Va) 
E E CfdMK*-,{aM)&^n{ah)Ci^i = ^^''^^MKic^h) (7b) 

E CfdM^%{ah)^t.{o^h) = E CfnK^'JAc^Pl) (7c) 

A more familiar form of this is the Kronecker relation of product reduction (8) ^ 0 ^-^^ 0 • • • . 
Another form is a body-to-lab coupling relation with M = fi + m and n = K — ft fixed in the /x or /2 
sums. The latter yields a sum over = K — n of body-fixed BOA waves Equation (5) giving lab-based 
$}(^ wave Equation (4). 
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(8) 



(9) 



(10a) 




(10b) 



Body-(un)coupling in Equation (10a) is an undoing of B OA-constriction by subtracting vi- 
bronic (/,/2) from {J,K) of BOA-wave ^j[u{e)] in Equation (10b) to make lab-fixed in 
Equation (10a) with sharp rotor quanta i? = J — I, J — I + 1 . . . or J + 1. In a lab-fixed wave 
of Equation (5) or (10a) rotor R is conserved but K and ft are not. A BOA wave $ j[i/(e)] of Equation (4) 
or (10b) has body-fixed vibronic K and that are conserved but rotor R is not. 

Note the following for Equations (8)-(10b). For Equation (8) we have constant (M = ju + m). Result 
Equation (9) is derived from Equations (6b) and (7c) with constant {n = K — ft). In Equation (10a) 
K = fl + n. In Equation (10b) M = ji + m. 

However, in both Equation (10a) and (10b) the internal bare-rotor body component 
n = K — ft is conserved due to a symmetric rotor's azimuthal isotropy. This n is a basic 
rovibronic-species quantum number invariant to all lab based perturbation or transition operators. Like 
a gyro in a suitcase, no amount of external kicking of the case will slow its spin. Only internal body 
operations can "brake" its n. 

The duality of lab v^*. body quantum state labels and external v^". internal operators is an important 
feature of molecular and nuclear physics, and it is to be respected if we hope to take full advantage of 
symmetry group algebra of eigensolutions. The duality is fundamental bra-&-ket relativity. For every 
group of symmetry operations such as a 3D rotation group R{3)iab = {• • • R(o^/37) • • • } there is a dual 
body group R{?>)hody = {• • • R(o^/37) • • • } having identical group structure but commuting with the 
lab group. Tensor multipole operators, discussed next, come in dual and inter-commuting sets as well. 
Generalized Duality is key to efficient symmetry analysis as shown beginning in Section 6.1. 

1.3. Mathematical Background for Tensor Algebra 

1.3.1. Unitary Multipole Functions and Operators 

Spherical harmonic functions y^((/)^) are well know orbital angular factors in atomic and molecular 
physics. They are special [n = 0)-cases of Wigner-L^^ functions Equation (3) as follows. 




where: [/] = 2/ + 1 



(11) 
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A diatomic or linear rotor must have zero body quanta (n = 0) and has a Yj^{(/)9) rotor wave. 

-matrix elements or expectation values of a multipole potential are proportional to Clebsch forms 
of Equation (7a). 



rf(0^oX,o ( W ^5 ( W ^;^*o 



(12) 



A multipole matrix is Equation (12) with factor depending on {f,k,J} but 

not {m^, m}. 



m' 



m 



(13a) 



Factor ( J'||v^|| J) is the reduced matrix element of and chosen by a somewhat arbitrary convention. 



[k] 



(13b) 



This particular choice simplifies bra-ket coupling and creation-destruction operator expressions for v^. 



'jk 

'mq 



J' 



m,m' 
=q—m 

= E^s^'.)( ' (-1) 

J J' 



J*\t 



J—m 



(13c) 



=q—m 



q m —m 



=q+m 



Other choices rescale eigenvalues but do not affect eigenvectors of a tensor or its transformation 
behavior Equation (14). (By Equations (7c) and (13c), transforms like Equation (6a) for a wave 
function (06^).) 

k 

vj = R(a/37)vjRt(a/37) = E ^l^'J^^^) (14) 

q=—k 

Examples of tensor matrices for J' = J = 1 to 3 are given in Table 1 . The J = 2 case is given in 
expanded form by Table 1. (Higher- J tables are g-folded to save space. Scalar (v^) = l/V[J]is left 
off each J-table in Table 2) 
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Table 1. Tabulated values for J = 1. 



(,:) 






V2 




f - ) 


1 


c ■) 


1 

2 








(■ ■ 

1 • • 

^. 1 .) 


1 

^/2 


/ 
\ 


'l . 
• 0 • 

/ • -ly 


1 




/. -1 • ^ 
• • -1 


1 










/I • A 
• 1 ■ 

V ■ V 


1 

^/3 









(v,^=-2...2)'= 
/l -1 l\ 
1-2 1 

V -1 V 


=1 

1 

1 


K=-l...l)"= 
/l -1 ■ ^ 

1 0 -1 

V 1 -1. 


=1 

\ • 

1 

V3 

1 A 




^1 ■ A 
• 1 ■ 

V- • V 


1 

^/3 



Historically, spinor J = 1/2 tensors shown in Table 3(a) are related to four Pauli spinor matrices 
and Hamilton quaternions {1, i, j, k} in Table 3(b) or Table 3(c). The latter appear in 1843 and are used 
for Stokes' polarization theory in 1867. The are U (2) algebraic basis of quantum theory for physics 
ranging from sub-kHz NMR to TeV hadrons and also applies to relativity. General U (k) algebra has fc^ 
generators Vq, v^, . . . , with a subset of k mutually commuting diagonal {q = 0) labeling operators Vq 
of the U (k) tensor algebras. The are related to elementary creation-destruction ejk = a] a/^ -operators 
and to their RES in the following sections. 
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Table 2. Unit tensor representations. 



K=-i...i>' 



1 1 

1 0 -1 

• 1 -1 



1 

1 
V2 



2 

72 



-V2 






1 


-^/3 




V3 


0 


-^/3 




%/3 


-1 






V2 



-V2 



1 

1 



3 -\/3 
2 

• V5 



1 



-V6 
V6 



-1 
V5 



-75 
-2 
V3 



-3 



1 

1 

\/28 



-1 1 
-2 1 
-1 1 



1 

1 

V2 

1 
x/6 



2 -x/6 72 
\/6 -1 -1 \/3 
\/2 1 -2 1 \/6 
-v^ 2 



\J=3 
2/ 



1 
1 

1 

Vl4 



5 
5 

75 



-5 
0 

/T5 

/lo 



75 

-7i5 710 

-3 -72 



72 
712 



-4 
-72 

/lo 



V12 

72 

-3 
-7l5 

75 



10 
l5 75 
5 
5 



0 



1 

V42 
1 

V84 
1 

^/84 



/ 3 \ J=2 . 

\%=-3...3/ 



1 


73 


1 


-1 




73 


-2 


72 


0 


-1 


1 


-72 


0 


72 


-1 


1 


0 


-72 


2 


-73 




1 


-1 


73 


-1 



1 

V2 

1 
V2 

1 

1 

^/To 







/^3 \J=3 








1 


-72 


72 -1 








72 


-1 


0 1 -72 








1 


1 


-10 1 


-1 


-1 


1 




72 


0 -1 1 


0 


-72 








72-1 0 


1 


-72 








1 -72 


72 


-1 


76 

1 












V6 



M=-4...4>^=' ^ 



\J=3 
4/ 



1 


-1 


73 


-1 


1 


1 


-4 


76 


-78 


1 


73 


-76 


6 


-76 


73 


1 


-78 


76 


-4 


1 


1 


-1 


73 


-1 


1 



1 

1 

V2 
1 

1 
1 

\/70 



3 


-730 


754 


-3 


73 






730 


-7 


732 


-73 


-72 


75 




754 


-732 


1 


/T5 


-740 


72 


73 


3 


-73 


-7l5 


6 


-715 


-73 


3 


73 


72 


-740 ^ 


/IE 


1 


-732 


754 


•75 


-72 


-73 ^ 


/32 


-7 


730 








73 


-3 


754 


-730 


3 








=-5. 


\J=3 _ 
.5/ ' 







1 

V2 
1 

1 

1 

V6 
1 

1 



1 

75 

1 

72 

1 
1 



-75 

-4 
-727 
-72 

-1 
0 

1 



1 

27 
5 

10 
0 

-1 
-1 



-72 
-72 

-7io 

0 
7To 

72 

72 



10 

5 

27 
1 



-1 

0 

1 
72 
727 

4 

75 



-1 
-1 

-V% 
-1 

-7^ 
-1 



1 

1 
1 
1 
1 
1 

a/84 



1 


-72 


1 


72 


-6 


730 


1 


-730 


15 


72 


-78 


10 


75 


-3 


715 


1 


-712 


3 


1 


-1 


75 



\ q= — 6. 
-72 

-78 
-10 
-20 
-10 

-78 

-72 



\J=3 _ 
6/ — 

75 
3 

715 

10 

15 
730 

1 



-1 

712 

-3 

•78 

730 

-6 

■72 



1 
1 

75 

72 

1 
1 



1 

1 

V2 
1 

1 

^22 
1 

1 

\/264 
1 

\/924 
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Table 3. Tabulated values and relation to quaternions, (a) Tabulated values for 
J = 1/2; (b) Simple Conversion from v to a; (c) Conventional quaternion- spinor relations. 



(a) Tabulated values for J=l/2 



-1 ■ 



-1 



-1 



-1 







1 




1 

V2 










1 

V2 



(b) Simple Conversion from v to <j 



-(J_ 



'VI 



(J. 



^+1 = +^+ 



1 



-1 



a,, = 



— z 



V2vg 
0 



(c) Conventional quaternion- spinor relations 




1.3.2. Tensor and Elementary Matrix Operators 



Coefficient 



J 



m 



J 



m 



of elementary operator 



m' / \m 



is the following CO or Wigner 3-j. 



J\ / J 



in I \m 



^\(^m'rn) where! q = m! — m (15a) 



m 



q m —m 



( ^\j'+j-k I [^] r<kJJ' 



qmm' 



(15b) 



Each matrix (v^)^ for J' = J = 1 to 5 is displayed in compressed form by the following tensor 
representation Table 2. 
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CG-3j relation Equation (13c) implies ^v^^ and (^em\m^ matrices have ortho-complete unit 
vectors of dimension d{J,q) = [J] — g = 2 J — g + 1 along g^^-diagonal of each [-/^]-by-[J] matrix. 
For example, quadrupole octopole and 2^-pole V2 share the g = 2 diagonal of J = 2 Table 2. 



^q=±2 



J=2 



\=±2 



^q=±2 



J=2 



2 / \J=^ 



• 1 



J=2 I 



e-2,0 



J=2 



• 1 




J=2 




3=2 



+\ -ie-i,) +A/-(eo,2 



+ 



+ 



+0(e_u 




J=2 




J=2 



• 1 



(16a) 



J=2 II / \ J=2 

eo,2 




3 / \ ^=2 




Tensor (v^\ relations easily invert to /e^/^^\ by inspection due to their being orthonormal sets. 



e-2,0 



eo,2 



J=2 



J=2 



J=2 



rj\^q=±2 



rj\^q=±2 



rj\^q=±2 



J=2 



J=2 



J=2 



+ M^{^q=±2 



/3/ 4 



'q=±2 



J=2 



J=2 



J=2 



J=2 



(16b) 



Any [J] -by- [J] matrix is a combination of elementary 



and thus also of (v^)^. This leads to 



RES maps that approximate [-7^] -by- [J] matrix (v^)^ eigensolutions by plotting related combinations of 

Yq^{e,(^) for select Oi^. 
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1.3.3. Fano-Racah Tensor Algebra 



Diagonal dipole- vector (rank k = 1) matrix y^l^ is seen in top row of Table 2 to be proportional 

to the angular momentum z-component matrix ^J^^ • Diagonal 2^-pole (rank-/c) tensors ^Vq^ are 
linearly related to powers = J^J^, = J^J^J^,... up to the fc^^ -power J^. This relates 

Vq ^ -eigenvalues to powers of ^ J^^eigenvalues m and, in turn, leads to an RES scheme to analyze 

J 



j eigensolutions . 

For example, matrix diagonals in Table 2 give elementary representations for J = 2. 



10(v^ 



(J=2) 



(J=2) 



(2) 



(2) 



11111 

2 10-12 



(17a) 



10 V, 



70 V, 



^0/ 



= [2 -1 -2 -1 2 
= (^1 -2 0 2 -1 
= f 1 -4 6 -4 1 



(17b) 



Powers of (jz^ in Equation (18) are combinations of ^v^^ found by dot products with vectors in 
Equations (17a) and (17b). 



11111 

2 10-1-2 
4 10 14 
8 10-1- 
16 1 0 1 16 



0\(2) 



^/5 



10 ,^.^(2) 



10 



10 



14 



(18) 



10 



70 



<v^> 



1\(2) 



<J°> 



■0\(2) 



, 4.(2) _ 3^70 / ox (2) 



1\(2) 



+ 



•2\(2) 



34V10 ^j,^(2) 



120 



(19) 



31^/70 ,j2x(2) 



(3) (7) (8) 



+ 24 ^-^^^ 



Triangle inversion of Equation (18) gives each (vq)^ in terms of powers (J^)^ = in 
Equation (19). RES plots depend on relating (vg)^ expansions Equation (19) in to Wigner (J, m) 
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[k][ in Equation (14) and Legendre polynomials ^qq{-6-) 

0 m —m I 



polynomials (—1)'^ 

Pfc(cos 6) in Equation (11) that are also polynomials of J 
shed light on the eigensolutions of the former. 



I J| cos 6. By plotting the latter we hope to 



2. Tensor Eigensolution and Legendre Function RE Surfaces 

Legendre polynomials occupy the central (00) -component of a Wigner-D-^ matrix. 

Examples of Legendre polynomials of cos^ = Jz/\J\ and = \J\ cos 6 are given below. 

1 



Pi (cos ^ 
Pa (cos ^ 

P3(cos^ 

P4(C0S^ 

Pc 

|J|^Pi(cos^ 

|J|2P2(C0S^ 
|J|^P3(C0S^ 
|J|^P4(C0S^ 



cos^ 



1 

2 



-cos^^ 
2 



cos 6 



cos^ e 



3 

8 

1 



30 



cos^^ 



+ ^cos4^ 

o 



J. 



u 

2 ^ 



2 ^ 



30, 



2 t2 



JrJ: 



(20) 



(21a) 



(21b) 



Classical Pt functions are compared with corresponding quantized (vq)"^ unit-tensor e- values in Table 4 
that generalize examples of tensor matrix ( J=l to 3) -eigenvalues in Table 2 and Equation (19) to any J 
and m = J, . . . , — J. The powers of m and J in (vg)'^, shown in Table 4 are taken to higher order in 
Table 5. 

Norm 2^^l^l\/2J + A; : -A; + 1 makes each (vg)'^ a unit vector. (Note: A ^ a : 



6 = (A + a)(A + a — 1) ... (A + 6).) In contrast, normalized Pt have Pa;(cosO) = 1. Coefficients 
of cos^ Q sum to 1 = Sc^. Square \c^\ usually do not sum to 1. 

Tensor values (vq)^, (vq)^, and (vq)^ in [. . . ]-braces of Table 4 equal Legendre functions Pq. A, and 
P2 in Equation (21b) exactly using J-expectation values Equations (22a) and (22b). However, for rank 
higher than k = 2, Pj. is only approximately equal to (vq)^ though the approximation improves with 
higher J. 



{Jz)t = m={\J\)icose, 

(|j,|) = y7(jTT)- J + J 



(22a) 
(22b) 
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Table 4. Forming (vq) from powers of J and m. 



( 



U / 7 




V / \/ [J] ^Omm 



2a/3 
V2J+2:0 

V2J+3:-l 

23v^ 
V2J+4:-2 



[1] 



For large- J values, the (vq )^ in Table 4 approach the P3, P4, . . . of Equation (21b) according to the 
relation: ^| ) [J{J + 1)]^^^. However, (vq)^ differ significantly from for low J. The 

classical Pk in Equation (21b) lack the small terms (—2/3, —5/6, etc.) that kill the (vq) in Table 4 
whenever J falls below strict quantum limits such as whenever J < \m\ or J < k/2. However, the 
quantum "killer" terms become negligible for larger J- values (J > k) and this makes tensor eigenvalues 
converge to Pk and thus to their RES plots. 

2.7. Angular Momentum Cones and RES Paths 

Quantum J-magnitude Equation (22b) introduces a quantum angular momentum cone geometry with 
quantized angles given by Equation (22a) as summarized here in Equation (23a) and (23b) for lab 
m = M and molecular body n = K. 



M 



cos 9m = , (23a) 



K 

cos el = , (23b) 

/7(7Ti) 

^ \ has sharp (zero-uncertainty) eigenvalue m orn on the lab or 

m, n / 

body frame z ox z axis, respectively. This sharp altitude and magnitude in Equation (22b) constrains 
vector J to base circles of cones making half- angle 6^;^ or 6^;[ with z ox z axes, respectively. Expected 
J -values appear in Figure 1 at intersections of quantized J -cones with the RES as explained below. 



An angular momentum eigenstate 



Int. J. Mol. Sci. 2013, 14 



Table 5. Forming (vq) from powers of J and m, expanded. 
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63 
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^(J + 3 : -2) 




^((J + 2 : 
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-1) 
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3(J + 2 : -1) 
f J(J+1) - f ) 


+ 




+ 2 : -1) 
5J(J + 1)+ ^) 
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35(J + 4 : -3) 
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RES energy level analysis begins by writing a multipole tensor expansion Equation (24a) of 
a general rigid rotor or asymmetric top Hamiltonian and then plotting the resulting surface using 
Equation (24a) 

H = AiJ.f + B{Jyf + C{J,f 

= ^-{A + B + C)T^+^-{2C-A-B)Ti (24a) 

+ i=(A-5) (T| + T!,) 



r° = (j,)2 + (j,r + (j,f = |J| 

= |JpP2(cOS^) 



3 2n 1 

- COS u 

2 2 



(24b) 



^ sin^ 9 cos 20 

Inertial constants (^4 = l/Z^, 5 = l/^y^ = l/h) combine J-tensor operators T^. Exact relation of 
(vq)^ and (vq)'^ in Table 4 to classical Pq and P2 in Equation (21b) is used in Equation (24b) for Tq 
andT^. 

A rigid spherical top {A = B = C) has only the term Equation (24a). Rigid prolate {A = B > C) 
or oblate {A = B < C) symmetric tops have only Tq and Tq terms with energy eigenvalues below. 



J 

^ m, n 



J 

m, n 

W 

o \m^n 



jjSymTop 

B{J,f + B{Jyf + C{J,f 



J 

m,n 
J 



(2B + C)r° + 2(C - B)T^ 



J 

.m, n 



B\J\' + {C - B){J,f 



m, n , 
J 

m, n , 

J 

m, n , 



(25) 



= BJ{ J+1) + {C - B)n^ 



Since a rigid symmetric-top involves only Tq^ and Tq, the ^;(-cones define its eigenvalues exactly by 
J-vector trajectories at angle-^;J where ^;J-cones intersect the following Tq^-RES shown in Figure 1. 



^^symTop(^) = ^(25 + C)T^{0) + ^(C - B)T^{e) 



(26a) 
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Inserting quantized-body cone relation Equation (23b) yields desired eigenvalues Equation (25) exactly. 

^^SymTop(^J) = ^(25 + C) J(J + 1) 

+ \{C-B)J{J+l)(lco,'ei-^^ 
3 V2 2y (26b) 

= J{J+l)\[{2B + C) + {C-B){K'-l)\ 
= J{J+l)B + {C-B)K^ 



Figure 1. J = 10 Symmetric top RES. Angular momentum cone of minimum uncertainty 
angle 0^^=17.55° intersects the highest K=J=1{) of the quantized J-path contour circles. 




Cone paths in Figure 1 are constant energy contours on symmetric top RES Equations (26a) and (26b) 
of constant J, They may be viewed as J-phase paths on which the J-vector may delocalize or "precess" 
on a circular O'^-cont around body z-axis. Or else one might view paths on Figure 1 as coordinate 
space tracks of the lab z-axis around the z-axis by assuming J lies fixed on the former. Either 
view describes J in the body-frame by Euler polar and azimuth angles —15, —7 with angle /3 = 9:^ 
and|J|2 = J(J+l)fixed. 

J = ( J^, Jy, Jz) = (— I J| cos7sin /5, | J| sin7sin /5, | J| cos /5) (27) 

The difference between quantum solutions and semi-classical Pk solutions can be easily plotted as in 
Figure 2. The figure shows a slice of the semi-classical surface and the uncertainty cones for each m 
from J to — J. The orange circles indicate the intersection of the surface with the uncertainty cones and 
the blue circles indicate the energy of the quantum value, (vq)^, placed along the uncertainty cone. 
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Figure 2(a) shows some divergence between quantum and semi-classical energies for (vq)^, while 
Figure 2(b) are exact for all J as in the J = 4 cases (vq)^ shown. As J and m are made larger 
than k then semi-classical values converge on exact eigenvalues as described below. 



2.1.1. Reduced Matrix and RES Scaling 

An RES is a radial plot along J direction —/3,—j that has hills where energy is high and valleys 
where it is low, but at all points the same magnitude | J| = ^/J{T^^. Origin-shift to keep RES radius 
positive and scaling to display hills, valleys, and saddles, may be needed to make useful RES plots. 

A scalar term s • Vq added to a tensor combination T = T^ktk^^Q does not affect the T-eigenv^c^c^r^" 
and neither does an overall scaling of T to cT. This is true since eigenvectors are invariant to adding a 
multiple si of unit matrix 1 to T or multiplying it by cl. (Of course, eigenvalues would, respectively, 
be shifted by s or scaled by c.) 



Figure 2. Quantum eigenvalues (blue) compared with semi-classical cone values (orange) 
for multipole tensor rank (a) = 4 (approximate), (vq)^ and ^4(008 B^) diverge for large 
m; and (h) k = 2 (exact), (vq)^ and P2(cos B^) correspond for all m. 






m=-3 




















(a) 



(b) 



Wigner-Racah tensor algebra defines a reduced matrix element ( J'||r^|| j) to serve as a scale factor 
for each Clebsch-Gordan tensor matrix element having Wigner-Eckart form Equation (13a). 



J' 



^ m 



rpk 

q 



J 



m J 



(28) 
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This matrix 



^0 



of quadrupole- J-tensor Tq 



Jq in Equation (24b) reveals some key points. 



2 ^ 2 













m / 


\m 



{cli^n.)-{J\\T'\\J) 



(29a) 



-m2--J(J+l) 
4[J] 

V2J + 3 : -1 V2 



rp2 



1 



m'--J(J + l) 



V2J + 3 : -1 



(29b) 



Reduced matrix element (J ||T^|| J) cancels norm factor 4-\/[J]/V2J + 3 : — 1 in Cq;^^^. The result 
is the quadratic Legendre form \J\'^P2{'m/\J\) found inside [. ..]-braces of Table 4 with norm 
4^/[k\/^/2J~+'3~!~^ outside the braces. (The latter is just a norm in Equation (29a) and (29b) 



multiplied by the factor y^^lTP] from definition Equation (15b).) 

Apparent conflicts in factors are due to having sum-of-i'^ware^Z-component normalization of unit v'^ 



on one hand and sum-of-component normalization of Pk on the other. Matrix elements 



J 



, m 



T 



or 




use the former since amplitude squares give probability. However, it is Mnsquared amplitude 



m / \m 

sum Ecfcthat measures anisotropy of a tensor T = 
(Each Pk contribues PkiO) = 1.) Expectation values 



Ecfc ■ Pk since Ec^ is a maximum T-amplitude. 
J 



, m 



T 



scale linearly, too, but J tensors may 



have extra scale factors. 

Tensor = in Equation (29a) and (29b) scales as | Jp = J(J + 1) and T'' 



J'^ scale as | J|*^. 

Factor 4:^/[J]/^/2JT3T^ of C^;^i^ reduces scale | to V2JTT = ^/[J]. Then the reduced factor 
(J ||r|| J) brings it back to |Jp. 



4[J] 



V2J + 3 : -1 



4(2J+1) 



^y{2J + 3)(2J + 2)(2J + 1)(2J)(2J - 1) 



(30) 



(J+l) (J + l)J(J-i) 



Each rank-/c part has a factor \J\^ = \J{J + 1)|^/^. Anisotropy of mixed-rank J-tensor T = • is 
T>\J\^Ck, and thus is quite sensitive to quantum number J. So also are the RES and related eigensolutions 
of T. 
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2.2. Asymmetric Top and Rank-2 RES 



Plotting RES of non-diagonal Hamiltonians for the asymmetric top Equation (24a) involves 2"*^ -rank 
tensors with reduced z-axial symmetry, nonzero g-values, and non-commuting Ja combinations. 
Each Ja in the quadratic expressions in Equation (24a) is replaced by its classical Euler-angle form 
in Equation (27). 

Or else, each tensor in Equation (24a) is replaced by a multipole fucntion = \J\''^q*{-, /S, 7). 
(Recall Equation (24b)). 



^^^AsymTop = ^ (| j| cos 7 sin /5) V 5 ( | J | sin 7 sin /5) V C ( | J | COS /5) ^ 



(31a) 



^Asy^Top ^ -^ + -^ + ^ ^0 + A B ^, ^ ^ + ^2^) 



^^^AsymTop ^ ^ + ^ + C p _^ 2C - A - ^ 



^0^ 



^-5 



76 



A + B 2C-A-B 2^ 

— 1 cos p H — sm p cos 27 



(31b) 



Forms Equations (31a) and (31b) of rank-(/c = 2)-RES are equal and give the same plots shown 
in Figure 3, but tensor form Equation (31b) reveals symmetry. Terms and (Figure 1) are 
z-symmetric and non-zero near z-axis while t^rms are asymmetric and vary as sin^ (5 with polar 
angle (5 between the J vector and the z-axis. As /5 approaches 7r/2, X^^ terms grow to give equatorial 
valleys and saddles in Figure 3 while X^ vanishes. 

Figure 3. Asymmetric top RES J = 10. 



Separatrix circle pair 
dihedral angle 




Asymmetric tensor operators T^^ are non-diagonal and do not commute with diagonal or with 
each other, and so i/AsymTop eigenstates as well as eigenvalues vary with coefficient {A — B) in 

Equation (31b). As n^i^^s symmetric-top states ' ^ 



into asymmetric-top eigenstates, 6^^ cone 



circles around the z-axis of Figure 1 warp into oval-pairs squeezed by nascent ovals emerging from the 
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X-axis and bound by a pair of separatrix circle-planes that meet at an angle 9^^^ on the y(orS)-axis. 
Figure 4 shows a range of RES and levels between symmetric-prolate top (B = A or 6^^^p = 0) and 
oblate top (B = C or 9^^^ = tt). A most-asymmetric case (B = C or 9^^^ = 7r/2) is midway between 
the symmetric cases. 

Figure 4. Asymmetric top energy levels with corresponding RES. 
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^^^P = tan- 



-1 l^-^l 
\B-C\ 



0 fox:B = A 
7r/2 for : 5 = (^ + C)/2 
TT for : 5 = C 



(32) 



As B first differs a little from A, off-diagonal ^"^^ asymmetric X^^ first "quench" degenerate 

J 



±i^-momentum eigenstate pairs 



into non-degenerate standing cos or sin-wave pairs. 



> = 



> = 




(33a) 
(33b) 



These states have nodes or anti-nodes standing on hills, saddles, or valleys of the RES topography at 
the principal body axes. Whether a wave is cos-like or sin-like at an axial point depends on whether it 
is symmetric or antisymmetric at the point and thus whether that point is an anti-node or node. Nodal 
location can determine whether a cos-like or sin-like wave has higher energy. 

As B differs more and more from A, off-diagonal T^2 will mix standing waves like |c^) with others 
such as |c^±2)' |^K±4)' Ic^ie) ^hat share the same i/AsymTop symmetry described below. 

2.3. Symmetry Labeling of Asymmetric Top Eigenstates 

Throughout the range of asymmetric cases in Figure 4 the symmetry of H^^^^ in Equation (31a) 
and (31b) is at least orthorhombic group D2 of 180° rotations R^^, R^, and R^ about inertial body axes 
that mutually commute (R^^R^ = R^ = RyR^), etc.). Unit square (R^ = 1, etc.) R-eigenvalues ±1 
label nodal symmetry (+1) or antisymmetry (—1) on each axis. D2 is an outer product of cyclic C2 
groups for two axes, say C2{x) and C2{y). x and y values also label nodal symmetry for the z axis 
since R^Ry = R^- A Cartesian 2-by-2 product of C2{x) and C2{y) symmetry character tables shown in 
Table 6 gives four sets of characters and four symmetry labels [Ai, Bi,A2, B2] for D2 = C2{x) ® C2{y). 

Labels {A,B) or (1,2) for (x) or (y) symmetric and anti- symmetric states follow ancient 
arcane conventions. We prefer a binary (O2, 12) notation for C2 characters and N-ary notation 
{On, In, "^n, • • • , {N — 1)n) for Cat characters S>^^{RF) where each label tun denotes "m-wave-quanta- 
modulo N" as in Table 7. 



^^^(R^) = e 



-im-p{27T/N) 



(34) 



This notation is used in correlation Table 8 between symmetry labels of D2 and its subgroups C2{x), 
C2{y), and C2{z), respectively. Each row belongs to a D2 species and indicates which C2 symmetry, 
even (O2) or odd (I2), correlates to it. The Table 8(a), 8(b) and 8(c) follow respectively from the columns 
R^, Ry, and of Table 6. An even (O2) D2 character is ^°^(R) = +1 and odd (h) is ^^""(R) = -1. 

J = 10 //A^y^Top.igygig p^gy^g 3 consist of two sets of five pairs [{Ai,Bi) (^2,^2) {Ai,Bi) 
(^2,^2) {Ai,B,)] and [(^2,^1) (^1,^2) (^2,^1) (^1,^2) (^2,^1)] separated by a single {A2) 
level. Each is related to RES a;-valley path pairs ^ [±10, ±9, ±8, ±7, ±6] or 2;-hill pairs ~ 
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[±6, ±7, ±8, ±9, ±10] separated by a single ^-path {A2 : Ky ^ 5), Even-K belongs to a (O2) column 
and odd-K belongs to a (I2) column of C2{x) Table 8(a) or C2{z) Table 8(c). 

Valley-pair sequence {Ai, Bi), B2) ... is consistent with (O2) and (I2) columns of the C2{x) 
Table 8(a), and hill-pair sequence (^2,^1), {Bi,A2) ... is consistent with (O2) and (I2) column of 
the C2{z) Table 8(c). This is because lower pairs correspond to x-axial RES loops of approximate 
momentum ^ ±10, ±9 • • • ± 6 while upper pairs correspond to z-axial RES loops of approximate 
momentum Kz ^ ±6, ±7 • • • ± 10 in Figure 3 

Table 8(b) for C2{y) is not used since ±^-axes are hyperbolic saddle points on one separatrix path, 
unlike the disconnected pairs of elliptic RES paths that encircle ±x-axes or ±z-axes at hill or valley 
points. Only a single E level exists in Figure 3 at the energy E'^^p of the saddle points and their separatrix. 

{j. =0 

^Sep ^ ^AsymTop^^^^ j^^ j^^ ^ ^j^j ^ ^^^^ J ^ |^|2 (35) 

[jz =0 



Table 6. Orthorhombic 4-group D2 = C2 x C2 character table construction. 











D2 = C2{x) X C2{y) 


1 




Ry 


R^ 


C2ix) 


1 


C2{y) 


1 R^ 


Ai = (0202).^ 


1 


1 


1 


1 


A = (O2). 


1 1 


X 1 = (02), 


1 1 


A2= (02l2)xy 


1 


-1 


1 


-1 


B = ih). 


1 -1 


2 = (I2). 


1 -1 


Bl = {l202)xy 


1 


1 


-1 


-1 


B2 — (l2l2)xy 


1 


-1 


-1 


1 



Table 7. Group character tables for cyclic groups of symmetry order N. (a) N = 2; 
(b) N = 3: e = e^'^/^. (j.) n = 4. 

(a) (b) (c) 



C2 


1 




C3 


1 


Ri 


R2 


C4 


1 


Ri 


R2 


R3 


(O2) 


1 


1 


(O3) 


1 


1 


1 


(O4) 


1 


1 


1 


1 


(I2) 


1 


-1 


(I3) 


1 


e* 


e 


(I4) 


1 


—i 


-1 


i 








(23) 


1 


e 


e* 


(24) 


1 


-1 


1 


-1 
















(3)4 


1 


i 


-1 


—i 
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Table 8. Symmetry correlation table between species of D2 and its axial subgroups, 
(a) C2{x) subgroup; (b) C^iy) subgroup; (c) C2{z) subgroup. 



(a) (b) (c) 



^2 ^ C2(X) 


(O2). (l2)x 


D2 D C2{y) 


(O2), (I2), 


D2 D C2{Z) 


(O2). 


(I2). 




1 


Ai 


1 


Ai 


1 






1 


A2 


1 


A2 




1 




1 


Bi 


1 


Bi 




1 


B2 


1 


B2 


1 


B2 


1 





turns into an i/AsymTop eigenstate expansion of states 




with K ±2p 



As symmetric H^y^ becomes a more asymmetric i/AsymTop pigm^e 4, a hill or valley path bends away 
from its ideal single-X symmetric-top cone circle at constant polar angle 0^ Equations (23a) and (23b). 

Each H^y"^ state ^ 

above and below K, and its RES path bends from constant 9'^ toward polar angles 6^x±2 5 ^k±4^ ^k±6 • • • 
above and below angle 9'^. At energy near the separatrix E'^^p, bending of hill and valley paths become 
more severe as they approach separatrix asymptotes where the polar angle range Equation (32) expands 
to 26^^^P or TT and the bend becomes a kink. 

It is conventional to label i/^^P eigenstate \E) by both and quantum values since \E) may use 
either a basis or else a basis. However, and do not commute. For energy E above E^^^, 
a \J,Kz) expansion is more compact and a dominate \Kz\ value may label \E). For E below E'^^p, a 
\J,Kx) expansion has a meaningful \Kx \ label. For E near E'^^p, X-labels are questionable. 



Table 9. Correlation tables between octahedral symmetric, O and various cyclic subgroups. 





(a) 








(b) 






0 D C4 


O4 I4 


24 


34 


ODC3 


O3 


I3 


23 


^1 ;C4 


1 ■ 






AUCs 


1 






A2iC4 




1 




A2iCs 


1 






E l- C4 


1 ■ 


1 




ElCs 




1 


1 


Ti l- C4 


1 1 




1 


TUCs 


1 


1 


1 


T2 -I C4 


1 


1 


1 


T2iCs 


1 


1 


1 



(c) 



(d) 



0 D C2(il) 


O2 


I2 


0 D C2(p.) 


02 


I2 




1 




^i;c2 


1 




^;c2 




1 


^2;c2 


1 




EIC2 


1 


1 




2 






1 


2 




1 


2 


T2iC2 


2 


1 


7^2 ;c2 


1 


2 
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Though a general form of the symmetry identification process may be unfamihar, it may implemented 
by computer. Group projectors Equation (36) can distinguish how each eigenvector splits with respect 
to subgroup operations. The product of these projectors and the calculated eigenvectors identifies the 
subgroup symmetry of each level. 



Only projectors in lower symmetry subgroups are used because they are easy to calculate and there 
are fewer in number. With the eigenvector projection lengths and knowledge of the correlation table 
between the molecular group itself and the subgroup one can start to deduce the eigenvector symmetry. 
As mentioned earlier, one correlation table is not enough to fully identify an eigenvector's symmetry, but 
using several subgroups one can assign symmetry. This process is simpler than calculating projectors 
of the full group, particularly if one can use a Cn subgroup and ^j"^ {g) in Equation (36) is found by 
Equation (34). 

This method can be significantly simpler than a traditional block diagonalization. Block 
diagonlalizing the Hamiltonian requires projectors of the entire molecular symmetry group rather than 
of the smaller subgroups. 

The disadvantage of this method is that it becomes unstable when clusters are tightest. As eigenvectors 
become more mixed with tighter clustering the algorithm may be unable to distinguish. Some RES paths 
and level curves indistinguishable to numeric projector then appear black. Symmetry definitions hold 
for asymmetric tops where J < 50. Spherical tops are quite challenging as seen in Section 7. 

2.4. Tunneling between RES -Path States 

N-atom inversion in ammonia, NH^, is an example of molecular tunneling modeled by a particle 
whose closely paired levels (inversion doublets) lie below the barrier of a double- well PES. An RES 
generalization, sketched in Figure 3, shows level pairs such as (Ai.Bi), (^2,^2), etc. as rotational 
analogs of inversion doublets. Here the tunneling between left and right positions on a PEs is replaced 
by an RES inversion between left-handed and right-handed rotation of an entire molecule. Instead of 
oscillation of expected position values (r) between PES valleys there is oscillation of momentum (J) 
between pairs of x-paths {+Kx ^ —Kx)in RES valleys or else between pairs of z-paths {+Kz ^ —K^) 
on RES hills. Section 7 describes this phenomenon in more detail for molecules of Oh and symmetry. 

3. Tensor Eigensolutions for Octahedral Molecules 

Section 2 has shown that asymmetric top molecules may be treated semi-classically, using only tensor 
operators and RES plots with a seperatrix between regions of local symmetry. Spherical top molecules 
experience such symmetry locality, but with greater variety of local symmetry. This section focuses on 
the added complication and convenience of higher symmetry as well as showing novel rotational level 
clustering patterns related to RES paths and tunneling. 




(36) 



9 
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3.1. Tensor Symmetry Considerations 



Theory of asymmetric top spectra in Section 2 may be generalized to a semi-classical treatment of 
tensor operators for or Oh symmetric molecules such as CH4 or SFq. The results contain level 
clusterings that first appeared in computer studies by Lea, Leask, and Wolf [20], Dorney and Watson [21] 
and followed by symmetry analyses [22,23] and others described below. 

Up to fourth order, any such molecule may be treated using the Hecht Hamiltonian [24] rewritten 
in terms of tensor operators below in Equation (37a) and (37b) that isolates the rank-4 tensor term 
Equation (37c). 




14 




12 



^0 




(37a) 
(37b) 
(37c) 



This is continued below to higher rank tensors with more comphcated structure [14]. The coefficients 
of each tensor operator may be found from a spherical harmonic addition-theorem expansion of points 
at vertices of an octahedron. Coefficients c„ „j are based on Equation (38), where is the spherical 
harmonic and f{r) is a position of octahehral vertices (100), (010), (00 — 1). 



JV 

A normalized sum of these coefficients gives the rank-6 Oh tensor as follows. 



(38) 



2^ ° 



(39) 



The first study of RES and eigenvalue spectrum with varying rank-4 and rank-6 tensor operators [25] 
expressed in Equation (40) revealed intricate level cluster crossing shown below. 



= cos(e)rW + sin(e)T[^l 



(40) 



Later studies [26] examined rank- 8 contributions such as expressed by Equation (41). Effects peculiar 
to combining Hamiltonian terms of rank- 8 and higher include extreme clusters. 



= cos((/))(cos(e)rW +sin(e)r[^l) + sin((/))r[^l 



(41) 



As with the asymmetric top Hamiltonian, the octahedral Hamiltonian uses non-axial operators shown 
in Equations (37b) and (39). Such operators involve more than Legendre functions, complicating 
purely semi-classical analysis. Thus, approximate solutions based on axial operators alone apply only 
asymptotically for high J > 10 and in regions away from RES seperatrices. Combined powers of J 
and Jz do not give all levels. Tensor operators provide a sparse banded Hamiltonian matrix, but full 
numerical diagonalization is needed to get all levels to high precision. 

Gulacsi and coworkers [26] explored how eigenvalues vary with and for J < 10 and small 
contributions of T^^l . Results below agree but extend to larger J and use RES topology. 



Int, J, Mol ScL 2013, 14 



740 



While the asymmetric top systems show clustering related to symmetry subduction from D2 to a 
related C2 subgroup, octahedral molecular clusters relate to a variety of subgroups. Oh may cluster into 
D4, Ds, D2, C4, Cs, C2 or other subgroups involving reflection or inversion. 

For simplicity, this discussion will focus on O rather than Oh molecules to make equations and 
correlation tables easier to display and interpret. Later Sections 6-7 give fuller discussions and explain 
the reciprocity relations that are behind correlation tables. 

For i^2- symmetric molecules, clustering patterns are described in terms of the correlation tables found 
in Table 8. The similar correlation tables for octahedral molecules are given for C4, C3 and C2 in 
Table 9. The columns of Table 9 represent the different clusters of rotational levels found within the 
spectra of given molecule at a given rotational transition. These clusterings are identified by their 
degeneracy as well as their RES location. Since symmetry labeling of octahedral group O differs form 
asymmetric top D2, a new coloring convention for O levels is defined: Ai is red, A2 is orange, E2 is 
green, Ti is dark blue and T2 is light blue. 

In the RES, rotationally induced deformation or symmetry breaking is seen from the shape of local 
regions of the RES involving a specific contour. Figure 5 shows two different RES plots, both globally 
octahedral, but with local regions that correspond to a subgroup symmetry of the octahedron. Figure 5(a) 
demonstrates a possible RES of an octahedral molecule with Hamiltonian parameters that allow for C4 
and local symmetry regions to be present. The C4 regions are identified by their location and by their 
square base. Similarly, the C3 regions are identified by their location and triangular base. In this case C3 
symmetric regions are concave while C4 regions are convex. This is not required and is dependent on 
Hamiltonian fitting terms that change the relative contributions of and T'^l Likewise, Figure 5(b) 
shows the C2 regions that are determined by their location and rectangular base. 

Cluster degeneracy is a hallmark of a specific symmetry breaking. While a symmetric top spectra may 
be resolved into rrij levels, a rotationally-induced symmetry-reduced spherical top has several identical 
z axes. The m j levels can then localize on a single symmetry-reduced local region. The number of these 
regions must equal the degeneracy of the cluster in that same region. This degeneracy, is also found 
using the sum of the numbers in the columns of Table 9 or by Equation (42) given °^ is the order of the 
molecular symmetry group and is the order of the subgroup. 



In the cases shown here cluster degeneracy becomes 6, 8 and 12 for C4, C3 and C2 respectively. 

3.2. Numerical Assignment of Symmetry Clusters 

As mentioned previously, it is possible to diagonalize the Hamiltonian and organize species by the 
order of each block, yet this alone will not distinguish all levels. For Hamiltonians defined by as 
Equation (37a) it is possible to analytically [25] determine the symmetry of each level. Once T'^l or 
terms are present, a numerical examination of eigenvectors is required to assign the symmetry of 
each level. Subgroup projectors are used here where the cluster degeneracy increases and the symmetry 
becomes challenging to distinguish. These projectors represent a simplification of the symmetry analysis 
of an octahedral molecule into projections onto C4 symmetric projectors. The correlation table for O D 
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C4, shown in Table 9, and Equation (36) give the information necessary for the assignment. Moreover, 
when using the subgroup C4 there are only four projectors to create and a clever choice of axis can force 
several of these projectors to be entirely real or entirely imaginary. Conveniently, the C4 projectors can 
be used to diagnose level symmetry for clusters in any subgroup region. 

3.3. Octahedral Clustering vs. RES Topography 
3.3.1. Variation of Topography 

The Hecht Hamiltonian Equation (37a) and its higher order analogues are generic Hamiltonians. Such 
Hamiltonians have numerous fitting constants specific to a given molecule and a given vibronic species. 
To better understand all such octahedral systems, one must focus on changes in the level spectrum and 
RES plots with varying contributions of T^, Tl^l and Tl^l 

7"[4,6,8] Equation (41) has two bounded parameters 6 and (f) so several plots are required to 
explore this parameter space. By setting T'^l contributions to zero the eigenvalue spectrum for T'^'^l 
in Equation (40) can be plotted for changing values of 9, relative contributions of 4^^ and 6^^ rank 
tensor terms. Figure 6 plots such an eigenvalue spectrum and also places the RES plots that go along 
with important parts of the level diagram and, conversely, points out what spots on the level diagram 
correspond to important changes in the RES plot. We note RES have circular ring separatrices not 
unlike those on D2 RES in Figure 3. 

To understand the behavior of the level diagram in Figure 6 it is critical to inspect the changing 
shape of the RES plots. In particular, the clustering of levels in the eigenvalue diagram is dependent 
on the localized symmetry regions of the RES at each value of 9. Locally, the RES forms hills and 
valleys of a lower symmetry than that of the molecule. The local symmetry must also be a sub-group 
of the molecular symmetry. Figure 5 identifies regions of local symmetry C4, C3 and C2 whose local 
rotation axis lie fixed normal to the RES at the center of each region, respectively, even as 9 varies from 
Figure 5(a) to Figure 5(b). For some 9 one or two of the Cn regions may shrink out of existence as shown 
below in Figure 6. 

Figure 5. Symmetry axes of RES for differing contributions of Tl^l and Tl^l. (a) 
and C4 local regions; (b) C2 local region. 




(a) 



(b) 
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3.3.2. Semi-Classical Outlines v^". Quantum Eigenvalues 

With this understanding of local subgroup regions it is possible to discuss more detail of Figure 6. 
The correspondence between the RES plots and the level diagram can also be seen by appending the 
eigenvalue spectrum in Figure 6 with the height of the C4, C3 and C2 axes. This serves two purposes: 
To confirm that the quantum spectrum sits inside the semi-classical boundaries and to see that there is a 
change in the eigenvalue spectrum corresponding to changes in RES topology. Figure 7 shows the same 
quantum spectrum as Figure 6, but also includes the height (energy) of each symmetry axis. The outlines 
are printed in bold and are labeled for which Cn axis they each belong. 

Figure 6. J=30 Energy levels and RES plots for T^^^^^vs.[4,6] mix-angle 9 with levels 
above 0=0° (extreme left), Tl^l levels at 0=90° (center), and -Tl^l levels at 0=180° (extreme 
right). C4 local symmetry and 6-fold level clusters dominate at 0=17° while C3 type 8-fold 
level clusters dominate at 0=132°. In between these extremes are C2 type 12-fold level 
clusters particularly around 0=80° where a C3 — C4 level-cluster-crossing of the top 14 
levels occurs. 
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Figure 7. Quantum spectrum of octahedral Hamiltonian (Equation (40)) with changing 9, 
Bold lines are the energy of the classical symmetry axis labeled. 




Section 2 described how to predict the error between a fully quantum mechanical calculation and 
a semi-classical approximation of the symmetric rotor rotational spectra. For the symmetric rotor this 
was done analytically. It is difficult to be as exact in calculating error for an octahedrally symmetric 
Hamiltonian, but a line plot can show when an RES plot fails to describe quantum mechanical behavior. 

Rather than plotting the Hamiltonian as Equation (40) we will arrange it as 

r[^'^i = (i-x)rw+xr[^i (43) 

This changes semi-classical outlines from cosines to lines and shows where quantum levels exceed 
semi-classical bounds and where an RES approximation fails. Also, x-line plots show by degree of 
avoided-crossing-curvature for each level the degree of its state mixing at x. 

The three plots in Figure 8 show these spectra and semi-classical outlines for J = 30, J = 10 and 
J = 4. Figure 8(a) shows that the quantum levels fit for all values of x at J = 30, while Figure 8(b) 
shows some small disagreement near x = 2 for J = 10. Figure 8(c) shows that for low J there is strong 
disagreement between quantum calculations and semi-classical approximations. 

Figure 8. Spectrum of Octahedral Rotor Showing Semi-Classical Boundaries Given 
Equation (43). (a) J = 30; (b) J = 10; (c) J = 4. 




(a) (b) (c) 



Int, J, Mol ScL 2013, 14 



744 



Indeed, such plots as Figure 6 have been created before, both for the RES plots and the level 
diagrams [25]. Next, we show such an analysis of T'^'^'^l and demonstrate how such a Hamiltonian 
can show a different type of topology than previously reported. 

3.3.3. Variation of T'^'^'^l Topography and Level Clusters 

The inclusion of eighth rank operators to the Hamiltonian dramatically changes the possible types of 
RES local symmetry and the related level clustering. While Figure 5 shows C4, C3 and C2 symmetric 
local structures for RES plots for T'^'^I Hamiltonians, Figure 9 shows a new kind of local Tl^'^'^l RES 
path pointed out there with Ci symmetry. (That means no rotational symmetry!) The path is repeated 24 
times and thus belongs to a single cluster of 24 levels. As shown in Section 6.7.2. the cluster spans an 
induced representation D^^ (Ci) t O, also known as a regular representation of O. 

Figure 9. RES with Ci local symmetry regions visible. 




Details of the two dimensional parameter space appear in a figure Table 10 containing RES 

plots for several {9, cp) points. To be consistent with Equation (41), the plots increase 9 from 0 to tt going 
left to right and cj) from 0 to tt going top to bottom. RES O levels are colored with the usual red for Ai, 
orange for A2, green for £"2, blue for Ti and cyan for T2. 

As expected from Equation (41), the top and bottom rows are opposites to one another. That is, where 
one RES has a hill (higher energy), the other has a crater (lower energy.) The RES at ^ = 0, 0 = 0 has 
convex C4 and concave Cs structure as does the RES at 0 = tt, 0 = tt, but opposite the shape of the 
RES at either 0 = 0, 0 = tt or ^ = tt, 0 = 0. The ordering of the levels is also opposite. These two 
extremal rows also have no eighth order contribution, so they produce simpler shapes than the others and 
are incapable of producing Ci local symmetry regions. The middle row shows a different behavior: all 
the diagrams are identical. Again, this follows from Equation (41) wherever 0 = 7r/2. 

While Table 10 shows only the RES plots along the parameter space defined by Equation (41), 
Figure 10 shows level diagrams with RES plots placed showing the symmetry and topology present 
at a given point in the space. The bold vertical lines next to the RES plots indicate the spot in the level 
diagram that particular RES plot would exist. Again, it is clear the 0 = 7r/2 case would be unchanged, 
so it is not shown. The 9 = n case is neglected as it is a mirror image of the 0 = 0 case. 
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Figure 10. Level diagrams of energy v^. 0 for given 0 with RES plots at selected positions, 
(a) 0 = 7r/4; (b) = 37r/4. 




(a) (b) 
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3,4, Criteria forCi Level Clustering 

Figures 9-1 1 show where the local regions of hills and valleys form on the RES depending on mixing 
angles cj) and 6, Unlike the local symmetry regions known previously, the local Ci structures associated 
with 24-fold level clusters have no rotation axis to locate their central maxima or minima on the RES. 
However, they do have bisecting reflection planes that must contain surface gradient vectors and an 
extreme point for which the gradient points radially. RES plots with Ci local symmetries are shown in 
parts of second, third and fourth rows of Table 10 as well as parts of Figure 10(a) and 10(b). Figure 11 
shows how Ci regions lie on hills or else valleys and how they can be arranged with their neighbors into 
either a square or triangular pattern. 

Ci clusters require tensors of rank- 8 with (j) between tt and zero as 9 varies. Momentum J must be 
large enough for its minimum-uncertainty ( J=X)-cone angle B j to fit in a Ci region. 

RES with J as low as J=4 may have Ci regions but fail to fit its (20|=56°)-wide cones. Ci clusters begin 
to appear around J=20 (26^^=26°) but even for J=30 (26^^=21°) are still barely formed in Figure 1 1(a). 
There a minimum uncertainty cone appears to barely fit within a separatrix on a Ci hill between C^ 
and C4 valleys of its ( J=30)-RES. Others are situated more comfortably in valleys of RES shown in 
Figure 1 1(b) where they appear to encircle a C4 axes as in Figure 10(a) where a corresponding cluster of 
24 eigenvalues in 10 levels appear at the lower left hand side of the level diagram. In Figure 11(c) they 
surround a C2 axis. 

With higher J and O/^-tensors of greater rank than A;=8, one expects clusters of 48- fold degeneracy 
corresponding to Ci -regions centered away from Oh symmetry axes or planes. So far, these are only 
beginning to be explored and analyzed. To analyze such complicated tunneling effects (and better 
understand older ones) requires an improved symmetry analysis developed in Sections 4-7. 

Figure 11. J = 30 RES for various rank /c=4, 6, 8 combinations giving Ci features, (a) Ci 
hills around C3 and C4; (b) Ci valleys near C4; (c) Ci valleys near C3, 
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4. Introducing Dual Symmetry Algebra for Tunneling and Superfine Structure 

For a system to have symmetry means two or more of its parts are the same or similar and therefore 
subject to resonance. This can make a system particularly sensitive to internal parameters and external 
perturbations and give rise to interesting and useful effects. However, resonances can make it more 
difficult to analyze and understand a system's eigensolutions. The tensor level cluster states give rise to 
spectral fine structure discussed in the preceding sections and that splits further into complex superfine 
structure due to J- tunneling that is the focus of the following sections. 

Fortunately, the presence of symmetry in a physical system allows algebraic or group theoretical 
analysis of quantum eigensolutions and their dynamics. Groups of operators {g,g' ,g" , ...) leave a 
Hamiltonian operator H invariant {g^Hg = H) if and only if each g commutes with it {gH = Hg). 
Then each g in the group shares a set of eigenfunctions with H. However, if (g^) and (g) do not commute 
then the {g') and (g) sets will differ. 

Hamiltonians may themselves be symmetry operators or linear expansions thereof. Multipole tensor 
expansions used heretofore are examples. Expanding H into operators with symmetry properties, such 
as (a'^a) or {T^), helps to analyze its eigensolutions since, in some sense, a symmetry algebra "knows" 
its spectral resolutions. The underlying isometry of a system's variables and states contains all the 
sub-algebras that are possible i7- symmetries. 

If i7-symmetry operators {g, g' , ...) also commute with each other [gg' = g'g, etc.) then all g share 
with H a single set of eigenvectors as discussed in Section 5. Such commutative or Abelian symmetry 
analysis is just a Fourier analysis where all H are linear expansion of its symmetry elements {g,g' ,g" , ...) 
and simultaneously diagonalized with H. Such g expansions define both Hamiltonians H and their states 
as described in Section 5. 

However, non-commuting (non- Abelian) symmetry operators {g,g' ,g" ,...) of H cannot both expand 
H and commute with H, This impasse is resolved in Section 6 by using a dual local operator group 
{g^g' iQ" 1 "') that mutually commutes with the original global group. Then local {g) expand any H that 
commutes with global g, while the global g define base states and their combinations define symmetry 
projected states. Roughly put, one labels location while the other labels tunneling to and fro. 

In Section 6, the dual group (^3 ^ C^^) of the smallest non-Abelian group {D^ ^ C^y) is defined. 
Dual symmetry-analysis is demonstrated for a trigonal tunneling system by group parametrization of all 
possible (L^a)- symmetric H matrices and all possible eigensolutions for each. The example shows how 
global {g) label states while the local {g) label tunneling paths. In this way symmetry labels processes 
as well as states. An added benefit is a kind of "slide-rule-lattice" to compute group products. 

In Section 7, the local symmetry expansion is applied to octahedral OcOh tensor superfine structure. 
Local symmetry conditions are used to relate tunneling paths to RES topography discussed previously 
and predict possible energy level patterns. The OcOh slide-rule-lattices appear in Figures 22-24. 

5. Abelian Symmetry Analysis 

An introductory analysis of tunneling symmetry begins with elementary cases involving homocyclic 
Cn symmetry of n-fold polygonal structure. But, it applies to all Abelian (mutually commuting) groups 
A since all A reduce to outer products Cp x Cq x • • • of cyclic groups of prime order. 
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5.7. Operator Expansion ofCn Symmetric Hamiltonian 

The analysis described here and in Section 6 deviates from standard procedure [27-31]. Instead of 
beginning with a given quantum Hamiltonian H-matrix, we start with a Cn symmetry matrix (r) and 
build all possible Cn symmetric (H) -matrices by combining n powers (r^) = (r)^ of (r) ranging from 



identity 



to inverse r 



n-l 



[32]. 



H = ror^ + rir^ + rsr^ + . . . + r^-2r''"^ + r^-ir^"^ 
= TqI + riY^ + r2r^ + . . . + r_2r~^ + r_ir~-^ 



(45) 



In Equation (45) the rotation r is by angle 27r/n so rotation is by angle n27r/n = 27r, that is, the 
identity operator = 1 = r^. Thus power-p indices label modulo-n or base-n algebras. If n=2, it is a 
Boolean algebra Ci C C2 of parity [+1,-1] or classical bits [0,1]. (U{2) spin-algebras of q-bits have 
47r identity but are not considered here.) 



Sum rule: p + p' = {p + p') mod (n) Product rule: p • p' = {p • p') mod 



n] 



(46) 



We construct the general H-matrix using Cn group-product tables shown below in a g"^g-form and a 
g^^g-form that is equivalent for unitary operators g^^ = g"^ . In each table the fc^^-row label g"^ matches 
/c^^-column label g so that the identity operator 1 = g~^g resides only on the diagonal. This example 
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(47) 
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The g^^g-form produces a regular representation i?(g) = (g) of each operator g as shown below. 
Each R{y^) is a zero-matrix with a 7 inserted wherever a appears in the g^^g-table. 
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(48) 



Int. J. Mol. Sci. 2013, 14 



749 



The Cn Hamiltonian (H) matrix has matrices from (48) inserted into expansion (45) of operator H. 
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(49) 



Matrices in Equation (49) are simply group tables Equation (47) with complex tunneling amplitude 
replacing operator r^. Parameters ro = (ro)* and = {r^Y match self-conjugate binary subgroups 
Ci C C2 = (l,r^) related by 1 = (r^)^. Both are real if matrix (H) is Hermitian self-conjugate 

{Hah = H^a)' 

Three distinct classes of tunneling or coupling parameters are depicted in Figure 12 using classical 
spring-mass analogs for quantum systems [22]. Tunneling matrices have a long history [33] going 
back to Wilson [34]. Here this is being revived to treat extreme J-tunneling and more recently by 
Ortigoso [17] and Hougen [35,36] to treat extremely floppy molecule dynamics. Both these tasks use 
tunneling parametrization that has so far been quite ad. hoc. To accomplish either of these tasks, or what 
will surely be needed, namely both tasks, we need a tighter symmetry analysis. The group operator 
scheme being introduced here seeks a way to achieve this. 

The P^-neighbor class has non-zero parameters ri=-r and conjugate r_i=-r*=-f coupling only 
nearest neighbors each with self-energy r^=Hi. The 2^'^-neighbor class has non-zero parameters 
r2=-s and conjugate r_2=-5*=-5 coupling only next-nearest neighbors with self-energy ro=H2. Finally, 
S^'^-neighbor coupling rs=-t=-t'' is real as required for binary self-conjugacy r^=(r^)^ 



Figure 12. Three classes of tunneling paths and parameters. 
(b) 2"^ Neighbor C. 
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5.2. Spectral Resolution ofCn Symmetry Operators 

Eigenvalues ^^ch operator are m^^ multiples of n^^ -roots of unity since all Cn symmetry 

operators g = satisfy = 1 and are characters of Cn symmetry operators. Magnetic or 
morfe-wavenumber indices m label a base-n algebra as do the power or position-pomt indices p in 
Equation (46). Spatial lattice points Xp = L-p{meters) are indexed by p while reciprocal- (A;) -wavevector 
space = 2Tim / L{per meter) is indexed by integer m. 

The Cn irreducible representations D^^\y^) as well as Cn characters. General group 

characters are traces (diagonal sums) of D-matrices (x^^Hs) = traceD^^\g)). Abelian group 
irreducible representations are 1 -dimensional due to their commutativity, and so for them characters and 
representations are identical. (x^^Hs) = D^^\g)) All this is generalized in subsequent Section 6. Any 
number of mutually commuting unitary matrices may be diagonalized by a single unitary transformation 
matrix. The characters in Equation (50) form a unitary transformation matrix ^ that diagonalizes each 
Cn matrix (r^). 

Tm,p = Xp/V^ (51) 

This T is a discrete (n-by-n) Fourier transformation. A 6-by-6 example that diagonalizes all matrices in 
Equations (48) and (49) and in Figure 12 is shown in Figure 13 by a character table of wave phasors based 
on i^(^)(r^) in Equation (50) or (51). The irreducible representations D^^\y^) or irreps play multiple 
roles. They are variously eigenvalues, eigenvectors, eigenf unctions, transformation components, and 
Fourier components of dispersion relations. This hyper-utility centers on their role as coefficients 
in spectral resolution of operators into idempotent projection operators p(^). p(^) are like 
irrep placeholders. 



n-l 



= E ^rP^'"^ = xlP^'^ + xlP^'^ + xlP^'^ + xlP^'^ + XJP^'^ + X^P^^^ (52) 
Equation (52) is column-p of Figure 13. Column-0 is a completeness or identity resolution relation. 

n-l 

^0 = ^ pM = 1 = pW + p(l) + p(2) + p(3) _^ p(4) _^ p(5) (^g3>| 



m=0 



Dirac notation for P^^^ is |(m))((m)|. Its representation in its own basis (eigenbasis) is simply a 
zero matrix with a single 1 at the (m, m) -diagonal component. P^^^ -product table in Equation (54) is 
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equivalent through Equation (52) to g-product table in Equation (47) but the P^"*) -table given below has 
an orthogonal (e.^.P^^^P^^) = 0) idempotent (e.^.P(^)p(^) = P(^)) form. 



ptp 

form 


p(0) 


p(l) 


p(2) 


p(3) 


p(4) 


p(5) 


p(0) 


p(0) 












p(l) 




p(l) 










p(2) 






p(2) 








p(3) 








p(3) 






p(4) 










p(4) 




p(5) 












p(5) 



\ 



(p(^)), 



(54) 



v / 

The location of each P^^^ in the P-table is a location of a 1 in its representation as indicated in the 
right hand side of Equation (54) in the same way that locations in g-table Equation (47) place I's in 
representations Equation (48). However, idempotent self-conjugacy (P^^ = P) makes row labels of 
P-table Equation (54) identical to its column labels, whereas only g = 1 and g = are self-conjugate 
in g-table Equation (47). 

Character arrays such as Figure 13 represent operator eigen-products between and p(^). 

r^pM _ ^^pM = pMr^ (55) 
Also character is the scalar product overlap of position state bra or ket with momentum ket or bra. 

Position bra: = {p\ = (0|r"^ 

Position ket: \xp) = \p) = r^|0) (56a) 
Momentum bra: {km\ = {{m)\ = (OlP^'^^Vn 

Momentum ket: \km) = \{m)) = P^'^^\0)^/n (56b) 
Momentum eigenwave ipkmi^p) character Equation (50) conjugated to e^^^^p and normalized by 

^kmM = {xp\km) = {p\{m)) = {x'^/Vriy 

= {{{m)\p)y = e^^-^^VV^ (57) 

Action of on m-ket |(m)) = is conjugate and inverse to action on coordinate bra {xq\ = {q\. 

^kmi^q-p- L) = {xq\rP\km) = (g|r^|(m)) 

= {q-p\{m)) = e-^^^^^(g|(m)) (58) 
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The same overlap results whether moves a (m)-wave ^^-points forward or moves the coordinate grid 
p-points backward. This Cn relativity-duality principle generalizes to non-Abelian symmetry and is key 
to operator labeling of coordinates, base states, Hamiltonians, and their eigensolutions. 

Figure 13. Cq characters (a) numerical table; (b) wave phasor table. 
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^P=0 
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-1 
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26 
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36 =-36 




-1 
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-1 
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-1 


46 =-26 








1 
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£2 


-1 


£2* 





Wavefunction: T'^(jc„) = 2'^ 




wave phasors 
1^ 2 3 4 



p(^) projects m-states with conjugate characters cp^ = (x^)* with factor 1/n so p(^)'s are 
idempotent and sum to 1. (SpP^^^ = 1) But, \km) has cj)^ with factor l/\/n to be orthonormal so 
its squares sum to 1. (Ep\{xp\km)\'^ = 1) Thus projection Equation (56b) of by P^^^ has a 
factor ^/n. Inverse spectral resolution Equation (52) sums over column points p using cj)^ from each 
row-(m) of Figure 13. Factor 1/n makes P^^^ complete (S^P^^^ = 1 in Equation (53)) and idempotent 
(p(m)p(m) _ p(m)) Equatiou (54)). 

P^"^ = C^"" j = (^"^' + + + + + Cr')/6 (59) 
First row ((m)=(0)-row) of Figure 13 is an average, i.e., sum of all symmetry operators weighted by 1/n. 



P(°) = I ) /n = (rV + + + + r^)/6 



(60) 



Thus factors ^/n = \/6 in state projections in Equation (56b) give state norms ^/n/n = l/^/n in 
Equation (57). 



pM|0) = ( ^ ) /y^ = (0-|O) + + 0-|2) + C|3) + d4) + C|5»/V6 (61) 

vP=0 
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The (O)-momentum or scalar state is a sum over the (m)=(0)-row of Figure 13 normalized by 

P(°^|0) = b) j /V^ = (|0) + |1) + |2) + |3) + |4) + |5))/y6 (62) 

5.3. Spectral Resolution ofCn Symmetric Hamiltonian 

Given Hamiltonian H expansion in Equation (45) in operators and the spectral resolution in 
Equation (52) of r^, there follows the desired spectral resolution of H. The eigenvalue coefficients 
^(^) of p(^) define the dispersion function cl;(/c^) of H in Figure 14(a) where it is conventional to 
center scalar origin (m)=(0). 

n—l n—1 n—1 n—1 n—1 

H = r^rf = J] rp J] XpF*^"^^ = J2 ^^""^P^""^ where : a;^"^) = ^pXp = ^(km) (63) 

p=0 p=0 m=0 m=0 p=0 

Positive fc^-axis Cq array [...(0), (1), (2), (3), (4), (5), ...] of Equation (54) shifts to a zone-center-array 
mod-6: [...(4), (5), (0), (1), (2), (3), ...]=[.. .(-2), (-1), (0), (1), (2), (3), ...] using Equation (46). 

Examples of dispersion relations for three classes of tunneling paths in Figure 12 are shown in 
Figure 14. Dispersion uj{km) for Cq symmetry depends sensitively on the Hamiltonian tunneling 
amplitudes for — 3 < p < 3 (or 0 < p < 6) in Equation (49), and for any set of eigenvalues 
u){km) there is a unique set of found by inverting Equation (63). 

n-1 

rp = J2 0>^'"V^ where : <j>^ = [x^T (64) 

m=0 

A common tunneling spectral model is the elementary Bloch P*-neighbor Sl(6)-model shown in 
Figure 14a, much like that developed in Reference [33]. For negative values of ri=-r, a Sl(6) spectra 
for Cq consist of six points on a single inverted cosine- wave curve centered at m=0 with its maxima at 
the Brillouin-band boundaries (m)=±3. This curve applies to Bl(n) spectra for Cn where n equally 
spaced (m) points lie on the dispersion curve between m=±n/2 for even-n. The n energy eigenvalues 
^(m) ^j,^ projections of an n-polygon. For n=6 that is the hexagon shown in Figure 14a projecting two 
doublet levels cj^^^^ and cj^^^^ between singlet cj^^^ and singlet cj*^^^ at lowest and highest hexagonal 
vertices as follows from Equation (63). 

^^'^"H^m) = roXo + rixT + r-iX-i = - 2r cos{2^m/Q) (65) 

The 2^^-neighbor 52-model (Figure 14b) has a two-cosine-wave dispersion curve. An equilateral 
triangle projects energy doublet levels [uj'^^\ uj^^^^ from its lowest vertex and a quartet [cj^^^^ cj^^^^J 
from its upper vertices. 

o^^^^''\km) = roXo + r2X2 + r-2X-2 = - 2s cos{47rm/6) (66) 
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Figure 14. Energy level dispersion for architypical tunneling parameters: Bl:ri = — r, 

B2:r2 = -s, B3:r3 = -t. 




The 3'"°' -neighbor i?3-model (Figure 14c) has a three-cosine-wave dispersion, which for n-Q and 
r2,=—t separates levels into an even-m triplet [uj^^\ijj'''^'^^\ below an odd-m triplet 

= roXo + rzx7 + r-zX-z = H,-2t (-1)'") (67) 

Combining of /c^^ -neighbor r/^ -terms gives dispersion cj^^^ as a fc-term Fourier cosine series that 
is, for real r/e, a sum of the preceding three Equations (65)-(67). However, real Vk imply symmetry 
that is higher than Cq, namely non-Abelian reflection-rotation symmetry such as Cq^ or D^^ and a 
corresponding degeneracy between cj^^^^ levels that will be treated shortly. Simple Cq symmetry allows 
six real parameters with complex ri and r2. Then Equation (63) implies six levels that are generally 
non-degenerate as shown in Figure 15. Complex ri = Irle^"^ of a ZBl model describes chiral magnetic 
or rotational effects that include Zeeman-like splitting of m-doublets. The projecting hexagon tilts by 
the "gauge" phase angle (j) = ti/12 as the ZB1(6) dispersion uj{m) shifts. Then m doublets (±1) and 
(±2) suffer splittings that are P^-order in (j) while singlets (0) and (3) undergo shifts that are 2^^-order 
in (j). 
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Figure 15. Zeeman shifted Bloch dispersion for complex parameter in ZB1(6) model: 
ri — re^^ with 0 = 7r/12. 




6. Non-Abelian Symmetry Analysis 

Characterization and spectral resolution in Equation (63) of a Hamiltonian H^^^^^ uses its expansion 
in Equation (45) in Abelian group Cq. Similar spectral resolution of a Hamiltonian H by a non-Abelian 
group G = [...gi, g2.-.] of non-commuting symmetry operators might seem impossible. To be symmetry 
operators of H, elements gi and g2 must commute with H, but that cannot be if H is a linear expansion 
of them like Equation (45). The impasse is broken by introducing operator relativity -duality detailed 
below. A Da-symmetric tunneling H with a 3-well potential sketched in Figure 16 is used as an example. 

6.1. Operator Expansion ofD^ Symmetric Hamiltonian 

The simplest non-Abelian group is the rotational symmetry = [1, r^, r^, ii, 12, is] of an equilateral 
triangle. is used to show how to generalize Cq operator analysis of the preceding section to any 
symmetry group. The analysis begins with a g^^g-form of group product table like Equation (47) for 
Cq. However, also requires a gg^^-form giving the same product rules but using inverse g^^ ordering 
|..r^, r^, ...|=|..r^^, r^^, ...| along the top instead of down the left side as is done for the g^^g-form of table. 
(The two ±120° rotations and are the only pair (r^^=r^) to be switched by conjugation). The three 
±180° rotations are each self-conjugate (ip^=ip) as is (always) the identity 1^^=1. 
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Over-bar notation is used for dual-group = [1 , r^, r^, ii, 12, 12] of "body"-based operators 
isomorphic to 'lab"-based group. 
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Figure 16. Rotation operators [1, r^, r^, ii, i2, is] for a symmetric square-well potential. 
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Matrix representations Equation (69a) for D'^ or matrices Equation (69b) for are given, 
respectively, by g^^g or gg^^-forms Equation (68) just as g^^g form Equation (47) for Cq gives matrices in 
Equation (48). 

(1) = (r^) = 

/i \ /. 1 . . . A 

1 



1 • 
• 1 



V 

(ii) = 



V 



1 • A 
• 1 
1 ■ 



/i 



1 • • 
• 1 • 
■ • 1 



\ 



V 



(ii) = 

/. . . 1 . A 

1 



1 • 
• 1 



(i2) = 



.1.7 



/. . . . 1 A 

1 

1 



1 



7 




(r^) = 

/. . 1 . . A 

1 • 
• 1 

1 



v- • 

(i3) = 



1 



• 1 7 



• • 1 

• 1 • 
1 • • 



1 



V 



r^) 

/. 1 . . . A 

1 

1 



(is) 



1 



1 7 



1 



1 • • 
• 1 • 



(69a) 



(69b) 



Int, J, Mol ScL 2013, 14 



757 



Most pairs of resulting matrices in Equation (69a) do not commute. (For example (r^)(ii)=(i3) 
does not equal (ii)(r^)=(i2).) Identical non-commutative product rules apply to the dual bar group 
matrices in Equation (69b). However, all matrices of the latter commute with all matrices of the 
former D^. This suggests that the Hamiltonian matrix, in order to commute with its symmetry group D^, 
is constructed by linear combination of bar group operators of [32]. 



H = rol + rir^ + r2r^ + iih + ^212 + ^313 (70) 

Ds symmetric (H) matrix Equation (71) generalizes Cq symmetric (H) matrix Equation (49) to a 
non-Abelian case. 
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(H) = 2^r^(g)= . . . (71) 

9=1 



6.2. Spectral Resolution ofD^ Symmetry Operators 

Spectral resolution of or any non-Abelian group G = [...gi, g2...] entails more than the Cq 
expansion into a unique combination of idempotent operators P'^=|a)(a| multiplied by eigenvalue 
D^'^\g) coefficients as in Equation (52). It is not possible to diagonalize two non-commuting gi 
and g2 in one basis since numbers (eigenvalues) always commute. If gi and g2 do not commute, 
their collective resolution must include tigtn-matrix coefficients D"^^^ involving nilpotent (N^ = 0) 
operators P^^n-\m){n\ well as idempotent (I^ = I) operators P^^rn-\m){m\ ^^^^ Equation (52). 

Unlike a commutative algebra of Cn idempotents, which are shown in Equation (54) and uniquely 
defined by Equation (59), a non-Abelian algebra yields a panopoly of equivalent choices of P operators 
that resolve it. The number and types of these P's is uniquely determined by size and structure of certain 
key commuting sub-algebras. The key to symmetry analysis of quantum physics is to first sort out the 
operators and algebras that commute from those that do not. It amounts to a kind of symmetry analysis 
of symmetry and leads to a far greater diversity than is found in commutative Abelian systems. 



6.2.1. Sorting Commuting Subalgebras: Rank and Commuting Observables 

The rank p{G) of a G-algebra is the maximal number of mutually-commuting operators available 
by linearly combining the operators g/e of symmetry group G. p{G) is also the greatest number of 
orthogonal idempotents P^ that can resolve the G-identity 1 as in Equation (53). {^G is total number 
or order of G. Here "^D^ and "^Cq both equal 6.) 

Cq rank is obviously equal to its order {p{Cq) = 6), but the rank of turns out to be only four 
(p{Ds) = 4). As shown below, can have no more than four P-operators that mutually commute 
though there exist quite different sets of them. On the other hand, has just three linearly independent 
P"-operators that commute with all of D^, and there is but one invariant set of them just as there is but 
one set of P^^^ for Cq in Equation (59). 
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Rank is a key quantum concept since it is the total number of commuting observables, the operators 
that label and define eigenstates. Of primary importance are G-invariant labeling operators that 
commute with all g and not just with other labeling operators. are uniquely defined within their 
group G and invariant to all g. (glG'g"^=lG) For example, total angular momentum and e-values 
J( J + 1) are i? (3) -invariant . 

Next in importance are labeling operators I/f^_2, belonging to nested subgroups of 

G=Hn in a subgroup chain GDi7^_iDi7^_2D...i^i. Multiple choices of chains exists since each 
subgroup link Hj. is not uniquely determined by the i//e+i that contains it, but each Ihj^ is invariant 
to all possible Hj^k at level- A; or below. 

For example, the z-axial momentum and its e-values belong to a 2^^ link in 
chain-R{3)DR{2z)DCQ{z). Given R{3) there are an infinity of R{2) subgroups besides the one for 
z-axis of quantization. or Jy are just two of an infinite number of possible alternatives to J^. Each 
i?(2^) has an infinite number of cyclic Cn{C) sub-subgroups. 

6.2.2. Sorting Commuting Subalgebras: Centrum and Class Invariants 



The centrum k{G) of a G-algebra is the number of all-commuting operators available by combining 
It is also the number of G-invariant P'^-operators. Students of group theory know /^(G) as the 
number of equivalence classes of group G. elements in Figure 16 are separated into three classes of 
elements [l],[r\r2], and [ii, is, h]- (i^{Ds) = 3) 

Elements in each class are related through transformation gi=gtS2S7^ by g^ in group G. Sum /^k of 
^c/e elements in g^^^'s class is invariant to g^ transformation. (It only permutes g^^^-terms in /^k thus /^k 
commutes with all g^ in G.) 



/^k where: /^k = J]] gj = ^Tsj, gtg/cgt 



(72) 



t=i 



The product table for L^a class algebra [/^i = 1 , = + , = ii + is + is] in Equation (73) 
below follows by inspecting group product tables in Figure 16 or Equation (68). It is a commutative 
algebra since each /^j commutes with each /^k as well as with each g^. This guarantees a class algebra 
has a unique and invariant spectral resolution. 
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(73) 



The first sum in Equation (72) is over the ^Ck elements in g/^'s class. (^Ck is order of /^k-) The second 
sum is over all ^G group elements. The number of elements gt that commute with gk ^^'^Sk, the order 
of gjt's self-symmetry Sk. Each group operator gk has a self-symmetry group consisting of (at least) the 
identity 1 and powers {gkY of itself. The order of class-fc is the (integer) fraction ^Ck=^G/^Sk. 
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6.2.3. Resolving All-commuting Class Subalgebra: Ccntmm=K.{Ds) = 3 

Spectral resolution gives class-sum operators /^i, /^2, and k^s as combinations of three i^3-invariant 
P^-operators with each of the /^k eigenvalues as coefficients. The /^3 characteristic equation found by 
Equation (73) gives three directly. 



0 = - 9/^3 = (/^3 - 3 • l)(/^3 + 3 • l)(/^3 + 0-1) 

= (/^3 - 3 • 1)P^^ = (/^3 + 3 • 1)P^^ = (/^3 - 0 • 1)P^ 
Standard notation ^1,^42, and E is used for the three invariant idempotents P"^. 



/^i = 1 • P^i + 1 • P^2 + 1 • P^ 
= 2 • P^i - 2 . P^2 _ 1 . 

/^3 = 3 • P^i - 3 • P^2 + 0 • P^ 



P^ 



( /^l + + /^3)/6 
( /^l + /^2 - /^3)/6 

(2/^1 - /^2)/3 



( 1 + + 
( 1 + + 



(74) 



ii + i2 + i3)/6 
ii - i2 - i3)/6 



(21 



r^)/3 



(75) 



Traces of matrices (gk) in Equation (69a) are zero excepting Trace{l) = 6. Traces of (P"^) then 
follow. 



traceP^' = 1, 



traceP^'^ 



traceP^ 4 



(76) 



This means (P^O and (P^^) ^j-^ q^q^ 1-by-l projectors while (P^) splits into two 2-by-2 projectors. 
The latter splitting is not uniquely defined until subgroup chain D^dCs or a particular DsDC2 chain is 
chosen, but relations in Equation (75) are invariant and unique. The hik coefficients inside parentheses 
of P"^ expansion give the character table for traces of irreducible representations (irreps). Irrep 
dimension £^ is trace of the a^^-irrep of identity gi = 1. 
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Xt = TraceD^ig,), r = Xi = TraceD^{l) 



(77) 



6.2.4. Resolving Maximal Mutually Commuting Subalgebra: rank = p{Ds) = 4 



Completing resolution of uses a product of two completeness relations, the resolution of class 
identity /^i = 1 in Equation (75) with the identity resolution of a subgroup C3 = [1, r^, r^] or else 
C2 = [1, 13]. In either case invariant P^ splits but P^^ and P^^ do not. In Equation (78) P^ is split by 
C2 into plane-polarizing projectors + P^ = Pg^^ + Pf^^^. 
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(l+ri+ r2+ ii + i2+ i3)/6 = P 
( 1 + + - ii - i2 - is)/ 6 = P^^pi^ 
(21 - - - ii - 12 + 2is)/6 = P^P^' 
(21 - - + ii + i2 - 2i3)/6 = P^P^^ 



(All other 

pQ;pm2 



(78) 
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In Equation (79) is split by Cs into Right and Left circular-polarized projectors 



T-\ f class algebra \ 
3 [completeness J 



^ I subgroup \ 
3 [completeness I 



Tj /group \ 
3 \completeness ) 

-I "p^i _|_ "p^2 _|_ "p£^ _|_ -vyE 

J- — -^OsOs ^OsOs ^Isls ^2323 



(All other 



-|_ _ p03 _^ pis _^ p23 

= P^3^03 = ( 1 + ri + r2+ii + i2+ i3)/6 = P^ip03 

= = ( 1 + + r^- ii- i2- i3)/6 = P^2p03 
where: ^ \ ^ o , ^1 P«P^3 = q 

Pf,L = Pf323 = ( 1 + + ^ )/3 = P^P'^ 



(79) 



In Equations (78) and (79), neither P^^ nor P^^ ^p^^ or change except to acquire some C2 or C3 
labels. The total number (four) of irreducible idempotents after either complete splitting is the same 
group rank noted before: p(i^3)=4. But, the i?L-circularly polarized pairs Pf and Pf^ split-out by 
C3=[l, r^, r^] differ from the linear xy-polarized pairs Pf^ and P^^ split-out by C2=[l, h]- ^x,x 
Vy y are, respectively, parallel (symmetric i3Pf =+Pf ) and anti-parallel (anti-symmetric i3P^=— P^ ) 
to x-axial 180^ rotation i3 in Figure 16 and will be used in examples. 

6.2.5. Final Resolutions of Non-Commuting Algebra: ^(£^3) = 6 

Mutually commuting algebras resolve into (I^ = I) operators ^Zi,m-\m){m\ ^hat sum to identity 
operator 1. They are split using the "one-equals-one-times-one" (1=1-1) trick in Equations (78) and (79). 

Non-commuting algebras resolve into idempotents and nilpotent (N^ = 0) operators P^^ri=lm)(nl 
that are split out using the following "operator-equals-one-times-operator-times-one" (g=l-g-l) trick. It 
is only necessary that 1 be resolved into rank-number p of irreducible idempotents as in Equation (78) 
or (79). (Here p{D^) = 4.) 

g = 1 . g . 1 = (P^;, + P^i + Pf;, + Pj;^) . g . (p;^, + + Pf,, + Pj;^) (80) 

The product in Equation (80) could have sixteen terms, but only six survive due to idempotent 

orthogonality Pj^jPf j. ^ ^^'^^j,k^tj' ^^^^ ^^^^ ^^^^ P^^ P^^ remain invariant and commute 

with all P'j j and all g. 



y->y ^ ' y,y o y^y 



This reduces to a non-Abelian spectral resolution of that generalizes resolution Equation (52) 
of Abelian Cq and includes two nilpotent projectors P^^. multiplied by off-diagonal irrep matrix 
components Dj^j^ as well as the four idempotents P^^ with their diagonal irrep matrix coefficients Dj'j 
that are not altogether unlike the L)(^)(r^)p(^) terms in Equation (52). ( Now X has matrix indices 

g= E Y.Y.DU9)nk (82a) 

irreps {a) j=l k=l 
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+D^Jg)P^^, + D^^^{g)F^^^ (82b) 
where :P« ■g-P^,= D« (g)P« P-^-g- P^, = Bf^,{g)P-, 

Terms {l/n)D^"'^*{rP)rP in Equation (59) of P^'") of C„ in Equation (52) generalize here to PJ^ and 
invert Equation (82a) to Equation (83). 

P-. = r/^G)5]i^-(^)g (83) 

Ds resolution in Equation (82b) has two irreps D^^ and D^^ of dimension ^^1=1=^^2 ^^id a third irrep 
of dimension £^=2 as noted in the first column of the character array in Equation (77). The irrep 
dimensions are related to the centrum tv(^Ds)=3, rank p(Ds)=4:, and order ^Ds=6. The following power 
sums of apply to any finite group G. 

G — centrum : tv{G) = ^ (i^)^ = Number of classes^ invariants^ or irrep types 

irrep{a) 

G — rank : p{G) = ^ {£^y = Number of mutually commuting observables ^g^^ 

irrep(a) 

G — order : ^(G) = ^ (^^)^ = Number of symmetry operators 

irrep(a) 

6.3. Spectral Resolution of Dual Groups and 

Spectral resolution shown in Equations (82a) and (83) of non-Abelian group G reduce g-h-product 
tables in Equation (68) to P-projector algebra. 

= <^"^4/PJ.. (85) 

Product tables in Equation (86) for projectors PJ^. generalize the Gq idempotent table in 
Equation (54). Non-commutativity entails a pair of tables like the g^^g form and gg^^ -forms in 
Equation (68) for "lab" g and "body" g operators. Tables in Equation (68) differ by switching conjugate 
pair and on side and top.(r^1^ = r^) The rest are self conjugate. (ii=ii, etc.) Similarly, tables 
in Equation (86) differ by switching conjugate nilpotent pair Pf^ and P^^. (PfJ=P^^) The rest are 
self-conjugate. (Pfj=PJj) 
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ptp 

form 


pAi 

XX 


pA2 

2/2/ 


p£;i p£;i 

XX xy 


TyEi TyEi 

^ yx yy 




XX 


pAi 

xx 














pA2 

2/2/ 








p^i 

p^i 
^ yx 






XX xy 

p^i p^i 

yx ^ yy 






p^i 

xy 

p^i 

2/2/ 








TyEi TyEi 

XX xy 

p^i p^i 
^ yx ^ yy 





ppt 

form 


pAi 

XX 


pA2 

yy 


p£;i TyEi 

XX yx 


p£;i p£;i 

xy yy 




pAi 

XX 


pAi 

XX 










pA2 

yy 




pA2 

2/2/ 








p^i 

XX 

p^l 

xy 






P^i 

XX 

p^l 

XX 


P^i 

xy 

p^i 
^ xy 




p^i 

yx 

p^i 

2/2/ 






p^l 

yx 

p^l 

yx 


P^i 

2/2/ 

p^i 

2/2/ 





(86) 



The p^^p and pp^^ tables in Equation (86) give commuting representations of projector PJ^ just as g^^g 
and gg^^ tables in Equation (68) give commuting (g)G-matrices in Equation (69a). Wherever PJ^. appears 
in a table, a "1" is put in its (p)-matrix. Putting "D^j^{gy' at each P^^ spot instead gives the following 
p^^p-representation (g) p of g since it is a sum of D'^kid)^^ Equation (82a). 



1 XX / 


\ yy / 




>Ei\ |p^i\ 
xy / \ yy / 

\ 


















\ 









(87) 



Conjugate pp^^ -representation (g)p of g has complex conjugate "Dj'^{gy' put at each PJ^ spot. The 
matrices in Equations (87) and (88) are transformations (g)p = T{g)GT^ and (g)p = T{g)GT^ of the 
respective matrices in Equations (69a) and (69b) by transformation T composed of D^j^{g) components. 
The Cq analogy is Fourier transform Equation (51) from Equation (48) to Equation (54). 
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(g)p 



Tt = 



I yx , 



I xy , 



I VV / 



V 



\ 



DlVig) 



Ding) 



D^ig) 



(88) 



Dl'Jig) J 



Matrices ...(f^)p, (ii)p, ... defined by Equation (88) commute with every ...(r^)p, (ii)p, ... defined by 
Equation (87) while each represents identical non-commutative product tables in Equation (68). 
Both use real [x, ?/]-based i3-diagonal irreps D^j^{g) given below. 



g 



ll 


l2 


l3 


1 


1 


1 


-1 


-1 


-1 


O 1 


/ ZL^ \/3\ 

1 O O 1 


( 1 • 



D^^ (g) 

1)^2 (g) 



D^x xy (g) 




(89) 



Appendix-A describes elementary derivation and visualization of D^j^{g) and their projectors P^j^{g). 



jk\ 



6.4. Spectral Resolution ofD^ Hamiltonian 

Hamiltonian H-matrix in Equation (71) has six parameters [ro, ri, r2, ii, ^2, ^3] or coefficients of its 
expansion Equation (70) in terms of intrinsic ^3 operators [1 = f°, f^, f^, ii, i2, is]. The parameters are 
indicated in Figure 17 by tunneling paths between the first base state |1) and other i^a -defined base 
states |g) = g|l) representing potential minima. 



Figure 17. Da-operator defined states and tunneling paths. 



r,= <l|H|l>=V 
= <r |H|l)=r/ 
r*= (r2|H|l)=7-2* 

= <ii|H|l>=i7 

i, = <i2|H|l>=^/ 

i, = <i3|H|l)=i/ 
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The resolution of H-matrix then follows that of g and (g)p-matrices. Any reduction of all 
(g)p-matrices, such as the [x, y]-reduction in Equation (88), also reduces the (H)p-matrix accordingly. 
Row-1 of (Hp) in Equation (71) has all six parameters. 



H. 



ah 



° G ° G 

^{m\g)D:;ig) = yr,D:;ig) 



9=1 



9=1 



(90) 



If the P-nilpotent pair are switched to ...Pfy, Pyx-' then (H)p and all (g)p (instead of all (g)p as in 
Equation (87)) are diagonal with eigenvalues H"^^ and i/^^ or block-diagonal with a pair of identical 
2-by-2 H^-h\ocks. 



p^2\ |p-Kl\ |P-^1\ |P-^l\ IP-^I 

yy / \ XX I \ xy I \ yx / \ yy 























(H)p = f (H)cft = 






XX 

Kx 


xy 

Ky 






\ ■ 






XX 

Kx 


xy 

Ky 



\ 



(91) 



The i/-block matrix components follow by combining Equation (89) with Equation (90). 



h: 



H. 



xy 



H. 



yy 



+ i3D^'*{i3)=ro+ri+rl+ii+i2+i3 

+ i-iD^''*{i3)=rQ+ri^-r{-ii-i2-i-i 
--roDllil) + nDllir^) + rlDHir'') + iiDlHh) + i2Dll{i2) + is^^xl^s) 
=(2ro-ri-rt-ii-«2+2i3)/2 

--roDllil) + TiDllir^) + r{Dll{T^) + ^lDJ(^l) + 12DH{i2) + ^s^^J (^s) 
-V?,{-n+rl-i,+i2)/2=H^: 

=roZ)^;(l) + riD^;(ri) + rlD^ir^) + ^lD,^;(^l) + i2Dll{i2) + ^s^^^Jl^s) 
=(2ro-ri-r*+ii+i2-2«3)/2 



(92) 



Irrep-dimension = 2 implies (at least) 2-fold degenerate E'-level since eigenvalues of identical 
i/^-blocks must also be identical, but only certain parameter values give diagonal i7^-blocks in 
Equation (92), i.e., real ri = r2 and equal ii =12- 



1 H^x 


xy 




\ Kx 


Ky ) 





2ro-ri-rl-ii-i2+2is V3{-ri+rl-ii+i2) 
V3{-rl+ri-ii+i2) 2ro-ri-rl+ii+i2-2is 
ro-ri-ii2+is 0 

0 ro-ri+iu-is 



(93) 



For: ri=r^ and: 11=112=^2 

These are the values that respect the local D C2 [1 , is] subgroup chain symmetry that gave (x, y)-plane 
polarized splitting in Equation (78). This is broken by a complex ri or by unequal ii and ^2- Complex 
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ri = \r\e'^'^ gives rise to complex rotating- wave eigenstates similar to ones in Figure 15 but, unlike that 
ZBl model, cannot split E'-degeneracy. Unequal ii and ^2 shift standing-wave nodes but cannot split 
^'-doublets either. £'-levels may split if H contains external or /a6-based operators g in addition to its 
internal or body-hastd g, but it thereby loses its symmetry. 

6.5. Global-Lab-Relative G versus Local-Body-Relative G Base State Definition 

Non-Abelian symmetry analysis in general, and the present example of D'^ resolution in particular, 
involves a dual-group relativity between an extrinsic or global "lab-based" group G=Ds on one hand, 
and an intrinsic or local "body-based" group G=Ds on the other hand. Each g in G commutes with 
each g in G. 

In the present example, the global "lab-based" group G=i?3=[l, r^, r^, ii, i2, is] labels equivalent 
locations in a potential or lab-based field and is a reference frame for an excitation wave or "body" 
occupying lab locations. 

On the other hand, the local "bod-based" group (5=53=[1, r^, r^, ii, i2, is] regards the excitation 
wave as a reference frame to define relative location of the potential or laboratory field. 

Quantum waves provide the most precise space-time reference frames that are possible in any situation 
due to the ultra- sensitive nature of wave interferometry. This is the case for optical coherent waves or 
electronic and nuclear matter waves. The latter derive their space-time symmetry properties from the 
former, and these are deep classical and quantum mechanical rules of engagement for currently accepted 
Hamiltonian quantum theory. 

Interference of two waves depends only on relative position as reflected in the following equivalent 
definitions of base kets for waves in a potential of Figure 16 with six localized wave bases 
[|1), |r^), |r^), |ii), |i2), \is)] in Figure 17. (We call this the "Mock-Mach Principle" of wave relativity.) 

|g^)=g^|l)=g,"'|l) (94) 

Key to this definition is the independence and mutual commutation of dual sets Equation (69a) 
and (69b). 

gjg/c = gkgj (95) 

Neither relation makes sense if we were to equate gk with g^^. The effect of gk is equal to that of g^^ 
only when acting on the origin-state The action of global '12 in Figure 18a is compared with local '12 
in Figure 18b that gives the same relative position of wave and wells. In Figure 18c product iii2=r has 
the same action as i2ii=r"^=r^ on 
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Figure 18. Ds-operators compared (a) Global i^, (b) Local i^, (c) 12 followed by ii. 




(c) i,i2|l)=i,|i2)=r|l)=r2|l) 

^\ / 



(Set 

' upright) 




Different points of view show how "body" g operations relate to the "lab" g. Starting from state |1), 



r^=r rotates lab potential clockwise (—120'') in a view where the body "stays put". The body wave ends 
up in the same well as it would if, instead, the body rotates counter-clockwise (+120'') by r=r^ in a lab 
frame that "stays put." 

In a lab view, effects of body operation g^^ and lab operation g^^ on |1) are the same except that g^^ 
also moves each body operation gj in the same way to g^gjg^^. The lab view of a lab operation g^^ does 
not see any of lab gj axes change location. The following generalization of lab-body relativity relation 
Equation (94) using Equation (95) shows how gj affects arbitrary \^k)' 



(96) 



= gfcgjil) = gfcgjg/git|l) = gitgjg/|gfc) 
6.6. Global versus Local Eigenstate Symmetry 

Applying projector PJ^ in Equation (83) to origin ket |1) gives a local-global symmetry-defined ket 

la \ 
\jk/' 



lo; 
\jk 



9=1 



(97) 
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The norm-factor N=^G/£^ is a non-Abelian generalization of the integral norm for Abelian Cn 
eigenket projection in Equation (61). Interestingly, the non-Abelian norm is also an integer since irrep 
dimension £^ is always a factor of its group's order ^G. 

A non-Abelian projection ket in Equation (97) has two independent symmetry labels j and k 
belonging to global-lab symmetry operators g and local-body operators g, respectively. Application 
of g-resolution Equation (82a) to ket Equation (97) is reduced by P-product rules in Equation (85) to 
the following global transformation. 



g|j,> = gpj,ii)^ 

= E E D^'k' (9) P'^'k'Pjk |1) ^ = E Df,^ (g) P^,, |1) ^/N 

j' = l k' = l j' = l 

j'=i 



(98) 



The corresponding local operator g first commutes through FJf, according to Equation (95) and is 
converted by Equation (94) to inverse global g~^ on the right of using Equation (82a) again. Finally, 
unitary irreps D°'{g~^) = D"^{g) are assumed. 



g \-,) = gpj, |i) Vn = pj,g |i) Vn = P|,g-i |i) Vn 

= E E D^'k' (9-') PJ.PJ^.^ |1) ^ = E D^k' (9-') |i) Vn 

j'=l k'=l j'=l 

= E D^,, {g-') \%,) = E Dtr, (g) \%,) 
j'=l f=l 



(99) 



A summary of the results is consistent with the block matrix forms in Equations (87) and (88). 

{%\ g |J,> = D^,^ (g) , g |J,> = DP, (g) (100) 

Choice of subgroup C2 = [1, 13] in Equation (78) leads to (x, ^)-polarized states (771)2 labeled by their 
eigenvalues (—1) 



f +lfor :j=x f +lfor :k=x (101) 

— Oj/ j < ^ — Ok' k ^ 

I — Ifor :j=y I — Ifor :k=y 

Physical significance of these global-(j') and local-(/c) values are now discussed using Figure 19. 

Wherever the global j is x or i3-symmetric (O2), then the entire wave is symmetric to x-axial rotation 
by TT in Figure 19a or horizontal reflection through the middle square- well in Figure 19b. Similarly, 
wherever the global j is y or ia-antisymmetric (I2), that is seen for each overall figure, too. 

However, if the local A; is x or ia-symmetric (O2), the local wave in each well has no node and is 
symmetric to its local axis of rotation by tt in Figure 19a or horizontal reflection of each square- well in 
Figure 19b. Similarly, wherever the local k is y or ia-antzsymmetric (I2), that antisymmetry and one 
node is seen in each well, too. 

Local and global symmetry clash along the ia-axis for states projected by nilpotent P^^ or P^^. The 
result is the x-axial wave nodes indicated by pairs of arrows in Figure 19. The \Eyx) wave in the lower 
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right of Figure 19b appears quite suppressed on the 13 -axis. However, the simulation of the \Exy) in the 
upper left seems to have its "node" coming unglued. 



Figure 19. L^a-symmetry waves (a) Sketch of projection; (b) 3-Well wave simulation 
(Compare with Figure 20). 




The "unglued" level is higher than and enjoys more tunneling. If tunneling increases so do 
parameters such as ri and r2 in Equation (92) that do not respect x-axial local subgroup C2 = [1 , is] . This 
breaks x-axial nodes and 13 local symmetry causing £'-modes to be less C2 -local and more like current- 
carrying above-barrier C3 -local waves rotating on r-paths. correlation arrays in Equation (102) with 
C2 or Cs indicate level cluster structure for extremes of each case. 



D3DC2 


O2 


I2 


1^3 ^Ca 


O3 I3 23 




Ai 


1 




Ar 


1 • • 




A2 




1 


A2 


1 • ■ 




E 


1 


1 


E 


• 1 1 





(102) 



Column O2 of array D C2 in Equation (102) correlates to Ai and E. The lower (Ai, £')-level 
cluster in Figure 19 has O2 local symmetry and lies below cluster-(742, E) that has local I2 symmetry 
according to the I2 column of Equation (102). Column O3 of table D C3 indicates that Ai and A2 
levels cluster under extreme Cs localization, but columns I3 and 23 indicate that each E doublet level is 
unclustered under Cs with no extra degeneracy beyond its own = 2). 

A classical analog of quantum waves states in Figure 19 is displayed in Figure 20 in the form 
of vibrational modes for an Xs molecule. A detailed description of this analogy in Appendix A 
includes modes of various local symmetry combinations analogous to those introduced above and in 
Sections 6.7.1 and 6.7.2 below. 
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Figure 20. 1^33(72(^3) -local symmetry modes of X3 molecule (Compare with Figure 19). 




6. 7. Symmetry Correlation and Frobenius Reciprocity 

The mathematical basis of correlation arrays in Equation (102) is a Frobenius reciprocity relation that 
exists between irreps of a group and its subgroups. This may be clarified by appealing to the physics of 
PJ^ -projected states such as are displayed in Figure 19 and by exploiting the duality between their 
local and global symmetry and subgroups. 

i^3- symmetric Hamiltonian H in (71) is made only of local g that couple states through local 
/c-indices by Equation (100) but leave all values of global j-indices unchanged. Thus a-eigenstates 
of H mix /c-values to form ^^-fold degenerate levels labeled by j-indices. (Recall = 2 equal 
sub-matrices Equation (93) in (91).) Further degeneracy or near-degeneracy ("clustering") occurs if 



Int, J, Mol ScL 2013, 14 



770 



inter-and-intra local tunneling coefficients decrease exponentially with quantum numbers thus isolating 
equivalent local modes into nearly degenerate sets of "spontaneously" broken local symmetry. 

In contrast to this clustering or "un- splitting" associated with local g symmetry operators, global g are 
associated with external or "applied" symmetry reduction that causes level splitting. Adding global g^ 
to a Hamiltonian H reduces its G-symmetry to a self-symmetry subgroup K=Sm consisting of operators 
that commute with g^. Adding a combination of g^ and g^ reduces K to an even smaller self-symmetry 
intersection group 5^ n 5^. 

Global g couple |J;.)-states through global j-indices according to Equation (98). The more global 
perturbations are added to a Hamiltonian H the more likely it is to split ^^-fold j-degeneracy (for £^ > 2) 
and/or linebreak alter eigenf unctions. 

6.7.1. Global "Applied" Symmetry Reduction, Subduction, and Level Splitting 

In the G=Ds example, adding matrix (r^) from Equation (69a) to (H) in Equation (71) reduces its 
symmetry to K=Cs=[l, r^, r^], and adding (ia) reduces it to K=C2=[1, h]- Adding a combination of (r^) 
and (is) completely reduces (H)-symmetry to intersection C3 n C2=Ci=[l], which corresponds to having 
no global symmetry. 

By reducing G to a subgroup KcG, each G-labeled a-level becomes relabeled by that subgroup K 
and split (if £^ > 2) in precisely the way that central G-idempotent is relabeled and/or split by unit 
resolution shown in Equation (78) or (79). The splitting in Equation (79) of idempotent into 
G3-labeled Pfgig plus Pf323 implies the doublet level uj^ splits into G3-labeled singlets uu^^ and cj^^. 
Both Ds singlets Ai and A2 end up relabeled with C3 scalar O3 labels. 

Ds D Cs P^relabel/ split D"^ relabel /reduce uj^ r el ahel / split 

A2 P^^ = P^^pOs = pA2^^ ^ JJA, I _ ^03 ^ ^A, _^ ^03 (103) 

^ _ p£;pi3 + p^p23 iCs^ ^uj^ ^ uj^^ 

Global D3DC2 relabeling and/or splitting is by Equation (78). Now Ds singlets have different labels O2 
and I2. 

-D3 D C2 P°'relabel / split D"" relabel /reduce uj'^relahel / split 

A2 P^2 = p^2pi2 = p^^2^^ ^ £)^2 I (72 ~ D^^ =^ OJ^^ OJ^^ (104) 

E = p£p02 + p£pi2 iC2^ ^ 

Center portions of splitting relations in Equations (103) and (104) use subduction symbols (|) to 
denote how each irrep-L^^ reduces to subgroup C3 or C2 irreps under their respective global symmetry 
breaking. Earlier studies [34] have referred to these multiple subgroup splittings as multiple frameworks. 
Each a-row of Equations (103) and (104) corresponds to the row a=Ai, A2, or E, of correlation array 
DsDCs or DsDC2, respectively, in Equation (102). 
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6.7.2. Local "Spontaneous" Symmetry Reduction, Induction, and Level Clustering 

Opposite to global GdK symmetry irrep subduction D'^{G)iK=...®d!'{K)®d!'{K)®... that 
predicts level- splitting is the reverse relation of local KdG symmetry irrep induction 
d!'{K)^G=...®D'^{G)®D^{G)®... that predicts "unsplitting" or level-clustering. In the former, an 
-dimensional irrep L)^(k) of global G-symmetry is reducible to smaller {l^ < block-diagonal 
irreps d^{k) of a subgroup K. In the latter, a K irrep is induced (actually projected) kaleidoscope-like 
onto coset bases of a larger induced representation d^'\'G of G that is generally reducible to G irreps D^. 

Base states \k tj ) of induced representation d^'l^G are each made by a G-projector PJ^ acting on local 
(i^-symmetry base state \k)=P^\k) defined by local K-projector P^. G-projection is simpler if PJ^ is 
also based on K-projection. (It helps to stick with one framework through this!) 

Of all L)3DG2 -projectors PJ^^.^ based on Equation (78), only P^^Os' ^^202' Pf202 ^^^^ ^^S^^ index 
/c2 = O2 . Only these can project induced states IO2 tjs) local base state IO2) corresponding to the 
02-column of D C2 array in Equation (102) having Ai and E. Similarly, A2 and E in the l2-column 
of Equation (102) correspond to Pf^i^, P^2i2' -^f2i2 projecting states II2 t%) from a local II2) 
state. Each projector P'j^k^ in Equation (104) has a G2-subgroup projector P^^ "right-guarding" the side 
facing each local ^2-ket {£2) = P^'^\h) that is similarly "guarded" by its own defining projector P^^ 
G2-subgroup projector orthogonality then allows only /c2=^2, giving the projection selection rules just 
described. 

Pi2fcJ^2> = Fl,^F''F'-\e2) = S'^'-F%eJ^2) = S'-'-\e2 "tl) (105) 

Each "right guard" projector P^ of PJ^ is part of a GdK subgroup splitting or subduction splitting 
D'^{G)iK=...®d^{K)®... as shown by DsiC2 examples in Equation (104). (These go back to the 
original DsDC2 subgroup chain resolution in Equation (78).) In Equation (105) each P^ selects which 
a-type induced bases \k tj ) and block-diagonal a-irreps can appear in a fc-induced representation 
d^{K)'\'G=...®D^{G)®..., and it implies a duality between induced (t) level-clustering and subduced 
(I) level-splitting as stated by the following Frobenius reciprocity relation. 

Number of D"^ in d^{K) t G = Number of d^ in i^^(G) ; K (106) 

The numbers on the left-hand side of Equation (106) would reside in the fc^^ -column of a 
GDX-correlation array such as in Equation (102) while the numbers on the right-hand side of Equation 
(106) would reside in the a^^-row of the same array. The examples in Equation (102) have only ones 
{1} and zeros {•}. A deeper correlation D^dCi to Gi symmetry, i.e., to no symmetry is a conflation of 
either the array DsDC2 or the array D^dCs in Equation (102) since Gi=G2nG3 is the intersection of G2 
and G3. 





Oi = li 




Ai 


1 




A2 


1 




E 


2 





(107) 



Gi local symmetry base |Oi)=|li) is the |1) in Figure 18 that contains scalar Ai, pseudo-scalar A2, 
and two E wave states in Figure 19 consistent with a single column of D3DG1 correlation array in 
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Equation (107). This column describes induced representation D^^ (Ci) t Ds, also known as a regular 
representation of D^. 

Reciprocity in Equation (106) also holds for non-Abelian subgroup irreps d^. is the smallest 
non-Abelian group so it has no such subgroups, but octahedral symmetry has non-Abelian and D4 
subgroups that figure in its splitting and clustering that are described in later Section 7. 

6.7.3. Coset Structure and Factored Eigensolutions 

Three pairs of kets in Figure 17 relate to left cosets [IC2 = (1, is), rC2 = r^C2 = (r^, ii)] 

one at each site. 

[(|l),|i3)), (|r^),|g)=rn|l),|i3)), (|r2),|ii))=r2(|l),|i3))] (108) 

Conjugate bras (g|=(l|g^^ relate to right cosets [C2=(l,i3), C2r^=(r^, 12), C2r=(r,ii)], again, one per 
C2-well site. 

[((lUisI), ((rl|,(i2|) = ((l|,(i3|)r^ ((r%(ii|) = ((l|,(i3|y] (109) 

C2 projectors P^'^=l{l-\-is)=P^ and P^2=2(l-i3)=P^ split bra (g| into ±-sum of bras mapped by left 
coset g^C2. 

P-^ =1 ((1| ± (isl), (r^l P'"^ =1 {{r'\ ± (isl), (r^l p-^ =1 {{r'\ ± (110) 

The same projectors split ket |g) into bases P^^ |g^ i^^^t are ±-sum of kets mapped by right coset C2g. 

[P-^ |1) =1 (|1) ± lis)), P'"^ \r') =1 iy) ± |i2)), P'"^ \r') =1 (|r2) ± (HI) 

S'-coefficients in i7-submatrix Equation (93) track C2 cosets. Row-(bra)-x 
terms in H^, line up in (+)-right-coset Ig+iag order . .(ri+ii), (r2+i2)- 
Row-(bra)-y terms in Hy, line up in (— )-right-coset Ig-iag order 
(ri-ii), (r2-i2)- Column-(ket) (±)-forms H^^ and H^^ line up in /e/t-coset order ...(ri±i2), (r2±ii). 
Either ordering gives the same matrix. Off-diagonal components H^y and i/^^ have x vs. y symmetry 
conflicts so coset parameters ± ^3) vanish. 

K]- K]y \ - [ + ^s)-l{ri+ii)-l{r2+i2) 0 • (ro + ^3) -f {ri+ii)+f {r2+i2) 



hL J V 0 • (ro - .3) +f (ri-zi)-f (r2-Z2) {ro-^s)-l{r,-^,)-l{r2-^2) 



<N K\v\ ) V 0 ■ + +f (ri+^2)-f (r2+.i) (ro-^3)-i(n-^2)-^(r2-^l) 

(112) 



y I \ ^ V U ^6) I 2 V -L ^iy 2 V' ^ V U ^cj; 2V' -L ^-L/ 2^ 



Kets P^|r^)=[P^|l), P^|r-^), P^|r^) span induced representation rf'^(C2)t^3, and P^|r^) span 
d^{C^'\D'^, Normalized states P^|r^)\/2 and P^|r^)\/2 correspond to a-type and 7r-type orbitals at 
vertex positions p=0, 1, or 2 in Figure 21. i^3 table in Equation (68) is reordered in Equation (113) 
below to display (^2(13) body-basis right-coset representation bra-defined by (g| = (l|g or ket-defined by 



Int. J. Mol. Sci. 2013, 14 



|l)=|g). The resulting i7-matrix in Equation (68) is Equation (71) reordered for cosets of C2 instead 
ofCs. 
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ro 


is 
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ro 


i2 
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r2 


is 


ri 


i2 


ro 





|ii) = 
iaf^l) 



13 

f2 



I2 
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(113) 



Figure 21. Induced representation C2t-D3 base wave states at vertex points p = 0,1, and 2. 
(a) 02tl>3 bases P^|rf)^; (b) lat^^s bases py\rP)V2. 




C2 ordered products in Equation (113) help reduce i7-matrix in Equation (71) to a direct sum of C2 
induced reps {d^^®d}^)'\D^ in Equation (114). Upper (02)-array in Equation (114) uses cr-orbital bases 
|r^) in Figure 21a while 7r-orbital bases |rp in Figure 21b span the (l2)-array. 
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- h 
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Any group component of Equation (114) or combination thereof is a possible tunneling matrix. 
Submatrices cf'^{g)'['Ds shown for g=r^, ii, and is reflect the effect of these operators on states in 
Figure 21a and similarly for (f'^{g)'\'Ds in Figure 21b. 



ri- 



• • 1 






1 ■ ■ 






■ 1 • 








■ ■ 1 






1 • • 






• 1 • 





(iiii) 



• • 1 






■ 1 • 






1 ■ ■ 








■ ■ -1 






• -1 • 






-1 • • 





«3- 



1 • • 






• ■ 1 






• 1 ■ 








-1 • ■ 






• • -1 






• -1 • 





(115) 



The O2 correlation in Equation (102) implies (f'^'l'Ds reduces further to irreps Ai®E that label 
the lower band of Figure 19. Meanwhile rf^^t^a reduces to irreps ^42®^ that label the upper band of 
Figure 19. Equation (91) shows ^10^2®^®^. 



7. Octahedral Symmetry Analysis 

Octahedral-cubic rotational symmetry O operations are modeled in Figure 22. Rotation inversion 
symmetry Oh=0 x Ci operations are modeled in Figure 23. In each case the larger g-symbols (such as f 1 
on top of Figure 22) label position ket states (such as |f i)=f 1 1 1)) while smaller g-symbols label axes of 
rotation in O (such as iQ on top facing edge of Figure 22 labeling that 180° rotation) or planes of reflection 
in Oh (such as the ax just above the z-axis on facing plane of Figure 23 labeling the x-plane reflection). 

Figure 22 is an "0-group shde-rule" since product ie • f 1 can be viewed as operator flipping a 
wave in position |fi) onto position |Rz), that is, i6|fi)=|Rz) giving product ie • fi=R^. Figure 23 is an 
"0/i-group slide-rule" (that does O products, too) and just as easily gives product • fi=S2 all without 
knowing what f 1 or §2 do. (As explained below, ri is 120° rotation about [111] axis and f 1 is its inverse 
located on the [III]-axis while f2 is on the [111] axis. §2 is f2 multiplied by inversion I-[lll]=[iiI]. ) 
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Figure 22. O operators distributed in cosets of C4 D C2. 




Note ie-transform of state |ri) (example: i6|ri)=|Ry)) differs from an ig-transform of operator ri 
(example: i6-ri-i6~^=r3^). The latter is divined easily by "slide-rule" as iq flips ri's axis onto ra^'s. 

Three Cartesian C4 axes of anti-clockwise 90° rotations R^^, R^, and R^ define directions [100], [010], 
and [001], respectively. Their inverses Rx=R^, Ry=R^, and Rz=R^ are also 90° rotations but around 
negative axes [lOO], [OlO], and [001]. A shorthand notation for 180° Cartesian rotations is pa;=R^, p^=R^, 
and Pz=^l' Trigonal Cs axes of anti-clockwise 120° rotations ri, r2, ra, and r4 lie along [111], [111], 
[111], and [111], respectively, while axes of inverses f i=r^, f2=r2, ^3=^3, and f4=r| lie along the opposite 
directions [111], [111], [111], and [111], respectively. 
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There are six C2 axes of 180° rotations ii, 12, is, 14, is, and ie located along [101], [101], [110], [110], 
[Oil], and [oil], respectively. This completes the five classes of O: [1], [ri..4, f 1..4], [pxyz]^ [^xyz, '^xyz]^ 
and [ii..6]. Including the rotations with inversion I yields five more classes of Oh- [I], [S1..4, §1. 4], [pxyz], 
[Sxyz^ ^xyz], and [ai..6] where Si..4=I • ri..4, [cr^yzM^ • Pxyz]^ [^xyz]=\^ • ^xyz], and [ai..6]=[I • ii..6]. c^'s 
are mirror-plane reflections in Figure 23. 

The "slide-rules" in Figures 22, 23 also help evaluate class products and construct left and right 
cosets of local symmetry subgroups. Three of the largest cyclic subgroups of O dirt tetragonal C4 
such as C4=[l, R^, R^=p^, R^=R^] displayed on the R^-face of the cube in Figure 22. In Figure 23 
the same face displays local symmetry C4^=[l, p^, R^, R^, 0-4, a^;, ^3, ay] that contains C4 plus pairs of 
diagonal mirror reflections [a4=I-i4, a3=I-i3] and Cartesian mirror reflections [ax=l-px, ay=l-py]. Each 
pair [ax, ay] and [a^, ^4] is a class as is rotation pair [R^, R^] or, singly, 1 and p^. The other five 
cube faces display cosets of the tetragonal subgroups C4yDC4 of O/^dO. 
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Figure 22 shows six 0-cosets of C4=[l, R^, p^, R^]. Opposite p^-f^ce has coset 

Px'Ca=[Px'>U, Py,h] in that order. The ri-face shows coset riC4=[ri, ii, r4, R^] in upper right of 
Figure 22, and the opposite r2-face has coset Y2'C/^=[y2, ^3, R^]. Top and bottom faces have cosets 
f r(:74=[f 1, Ro:, fa, ie] and YrC^=[Y2, R^, y^, is]. 

Each g-C4-coset element g-Rf {p = 0..3) transforms the 1-face to the same g-face and orients it 
according to a C4 element R^ as it permutes the list of its elements accordingly. Each face may be 
labeled by any element g-Rf in its coset. An i-class labeling by 1, i3(or 14), ii, 12, ie, and is of C4 cosets 
in Figure 22 is as good as any other. 

Figure 23 shows six O/^ -cosets of C/^y (counting C/^y itself) in a geometric display that also shows 
eight trigonal cosets of C^yZ:>C ^-[111] and twelve dihedral cosets of C2x;DC2-[101]. Figure 24 shows 
three symmetry points of Figure 23 forming a triangular cell with sides that are on reflection planes. 



An order-8 axial symmetry C^^ lies on the tetragonal-z- [001] -axis of a cube face or octahedral vertex. 
An order-6 C^^ lies on the trigonal- [111] -axis of a cube vertex or octahedral face. Finally, there is 
a dihedral-C2^ [110] -axis of a cube or octahedral edge. Lines between the axes have bilateral local 
reflection symmetry Cy{y)=[l,(jy], C^(2)=[l, ^2], or C^(4)=[l, 0-4], fundamental symmetry operations 
whose products generate all others. Figure 24 is like a reduced Brillouin Zone of the Oh lattice. 

Each subgroup spawns a coset space and a set of induced representations of full Oh symmetry that 
generalize the C^^ induced representations in Equation (115) and base kets sketched in Figure 21. 
Correlation tables between O or Oh and its subgroups LdG tell which O or Oh irreps, states, and energy 
levels arise from each coset space. As local symmetry reduces and its order °L decreases, the coset 
dimension d=°GrL grows proportionally with a corresponding increase in number of irreps and levels 
in LtG-induced representation cluster spaces. Examples are given below for G=0 and in Section 8 for 



Figure 24. Oh local symmetry (a) C^^; (b) Cs^; (c) C2^. 




3v 



G=Oh. 
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7.7. Octahedral Characters and Subgroup Correlations 

Spectral class resolution of O generalizes that of D3 in Equation (75) to give character array 
Equation (116). 



0 group 


g = i 


ri_4 

fi_4 


Pxyz 


Rxyz 
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-1 
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Ti 


3 
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-1 
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0 


-1 


-1 


1 



Cyclic subgroup C4(R^), C'^{y\), and C2 characters correlate to O according to arrays in Equation (1 17). 
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T2iC2 
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(117) 



Equivalent subgroup correlations 0Di7 and ODgi/g"^ share elements in the same 0-classes and have 
one correlation array. Thus all three C4 local symmetries have one correlation table in Equation (117), 
as do all four C3 subgroups. However, Oz^C2{pz) and 0DC2(ii) correlations differ since ii and pz have 
different 0-class and characters in Equation (116). 

Projectors PJ^ and irreps D^j^ of O depend on choice of local symmetry just as L^a projector splitting 
in Equation (78) or (79) depends on choice of correlation DsDC2 in Equation (104) or D^dCs in 
Equation (103), respectively. Sub-labels (j, k) range over C2 values [O2, 12] or else C3 values [O3, 13, 23] 
while a tetragonal correlation ODC4 will use sub-labels (j, /c)= [O4, 14, 24, 84]. 

The m4 or else ms unambiguously defines all O states since no ODC4 or ODC3 correlation numbers 
in Equation (117) exceed unity. However, 0DC2{h) correlations cannot distinguish all three sub-levels 
of Ti or T2 wherever a number 2 appears, and the OdC2{Pz) correlation leaves the E sub-levels 
unresolved, as well. A full Oh labeling resolves the first ambiguity as shown below, but we consider 
the unambiguous ODC4 case first. (C4 resolves C2{pz) ambiguities.) 
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7.1.1. Resolving Commuting OdC^ Local Symmetry Subalgebra: Rank = p{0) = 10 

The C4 correlation table in Equation (1 17) shows how invariant class projectors (expanded below 
in terms of O characters x^g tablt shown in Equation (116)) will split into irrep projectors P^^^^ 
when hit by C4 local symmetry projectors p^^. The latter are expanded in terms of C4 operators 
weighted by character eigenvalues = {x^^Y using Equations (57) and (59). 
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pa 


(P04 


+P14 


+P24 


+P34) ■ 


1 • 




-^0404 
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P^ = 
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_^pE 
^-^2424 
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1 pT2 

^-^2424 


+0 

,pTi 

_|_pT2 

^■■^^3434 



(118) 



The five class projectors P^ are 0-invariant and commute with all twenty-four 0-operators 
(1, ri, r2, ...is, ie). So do the five class operators (/^o, ^r^^ ^p^^ '^Rk^ '^ik) which each P^ is expanded as 
follows. (Recall D'^ classes in Equation (75).) 



^ 2^ Xt^k = where: a = A^, ^, E, Ti, or 

k=0 

24 



[xo 1 + + r2 + .... + f4) + + + Pz) 



+X^^(R, + R, + .... + R,) + xliSi + h + .... + ie) 



(119) 



Each of the irrep projectors P^^^^ is obtained from its invariant P^ by product P^p^4=p^4P^ 
following Equation (118) with each of four C4 local symmetry projector p^^. 



^ p27rim-p/4 
p=0 



P04 =(l + R, + p, + Rj/4 
Pi, = (1 + iR, -p,- iRj/4 

P24 =(l-R, + P,-RJ/4 
P3, = (1 - «R, - + iRj/4 



(120) 



As the five (O-centrum-5) projectors split into ten (O-rank-lO) sub-projectors P"^^^, the five O 
class sums Kg split into ten Ci-invariant sub-class sums Cfc(fc=1..10). 



°0 



k=0 

where: = (121) 

The resulting ten products -^Pn^n^ listed in Equation (122) of diagonal 

irrep coefficients D^^^^{gk) in terms of twenty-four group elements gk that have 
been sorted into ten sub-classes that have C4{z) local symmetry. The ten 

irrep projectors Pn^n^ ^4 local-invariant, that is, they commute with four 

C4-operators (1,R^,R^ = Pz,R^ = Rz) but not the whole O group like the P^ do. The ten 
sub-class-sum operators Ck, into which each P^^^^ is expanded in Equation (122), are each individually 
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invariant to Rf, that is RfCjt=c/eR^, and D^^^^{gk) is the same for all gk in sub-class c/^. Note that a 
sum of rows belonging to P^^^^ between horizontal lines in Equation (122) yields corresponding 
character values Xk-^^^^^^^idk) 0-character array Equation (116) and effectively "unsplits" the 
sub-classes. 
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+5 0 
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(122) 



Without evaluating Equation (122), one may find ten sub-classes by simply inspecting 

Figure 22 for operations in each 0-class that transform into each other by C4 operations Rf only. The 
0-class of eight 120° rotations r^^ split into two sub-classes, one [ri, r2, fa, f4] whose axes intersect four 
corners of the +z front square, and the other [fi,f2,r^,r^'\ whose axes similarly frame the —z back 
square. The class of six diagonal 180° rotations ik split into a sub-class [11,12, ^5, ^e] whose two-sided 
axes bisect edges of the Iz squares, and sub-class [^3,^4] whose axes are perpendicular to z-axis and 
bisect edges of Ixy side squares. The 180° rotational class [p^, py, p^] splits similarly into sub-classes 
[px, Py] and [pz] with axes perpendicular and along, respectively, the R^ axis. The 90° class splits, as 
indicated in the top row of Equation (122), into a sub-class of four perpendicular xy-axml rotations and 
separate sub-classes for and Rz. 

The inverse to Equation (121) expresses the ten subclasses in terms of the ten diagonal irrep projectors 
using the same (albeit, conjugated) array of D^^^^{gk). However, column and row labels must switch 
and acquire different coefficients. 

7.1.2. Resolving D-matrices with C4 Local Symmetry 

Off-diagonal D'^^^^^gj.) matrices derive from products of diagonal irrep projectors in Equation (122) 
using Equation (82b) repeated here. 



(124) 
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Scalar Ai and pseudo-scalar A2 are given first then E, Ti, and T2 irrep matrices for the fundamental 
ifc-class of O. 

Do,'o,{ikik) =ii + i2 + i3 + U + i5 + i6 ^^25) 
-^•^424(^^11^) = -{ii + 12 + is + ii + ib + ie) 



D'^iikk) = O4 
24 



+ 12 + i5 + k) -\- k + H 
h.c. 



"fiii + i2 



«5 



5(^1 + «2 + ^5 + ie) 



16) 
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(126) 
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I4 




34 








24 




I4 


+ ^{k + k + k + k) 


mk 


+ k- 


k 


- k) 




-k) + ^{k 


- h) 








-i{k - 


k) 










34 


h.c. 


+Uk 


+ k + 


k 


+ k) 




-k) + ^ik 


- k) 


O4 


h.c. 




h.c. 








+{k + k) 





Symmetry of C4CO subclass [^1,^2, ^5, ^e] and [^3, i^] would demand equality of parameters for each. 

anrf, ^3 = ^4 = ^34 = ill (128) 



^l = ^2 = ^5 = ^6 = H256 = ^i, 



Setting each parameter to the inverse of its sub-class order (ik=l / )) reduces each matrix to diagonal 
form and gives the diagonal D^^^^{gk) given in Equation (122). Classes r, p, R behave similarly. 

7.1.3. Resolving Hamiltonians with C4 Local Symmetry 

An octahedral Hamiltonian H = Ylk=i 9kSk with local C4{z) symmetry is resolved by sorting gk into 
its C4{z) sub-classes Ck and then into P^^^^ whose coefficients are the desired H eigenvalues e^^. Zero 
off-diagonal H^^^^ = 0 and C/^-local symmetry conditions shown in Equation (128) arise from 
Equation (122) consistent with Figure 22. Tunneling parameter H256=^i 

from +z-axis to its P^-neighbor ±x or ±.y axes may dominate flip-tunneling 
^34 = ill to 2^"^ neighbor-z-axis. The i-columns of Equation (122) 

(or matrix diagonals in Equations (125)-(127)) give ii and in contributions to eigenvalues 
fisted in the i^-column of Table 11. Clusters {e^l , eg , e^J through ( 

^34 5 ^34 ) plotted in Figure 25 for 

select values of parameters ii = ii256 and in = ^34. 

One expects the parameter in for 2^^-neighbor tunneling to be exponentially smaller than ii for 
adjacent tunneling so the {in = 0)-cases are drawn first in Figure 25. While the z-class operations 
are most fundamental (all operations are generated by products of ik) other operations also generate 
P^-neighbor transformation. Three class parameters i?^y(90°), ri(120°) , and ii(180°) label P^-neighbor 
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inter-C4 axial tunneling paths that have the same zj-level patterns and splitting ratios as (zn=0)-cases in 
Figure 25 but with differing sign. (Signs differ since each sub-class eigenvalue set must be orthogonal to 
all others as shown below.) Level patterns in Figure 25 are reflected in spectral patterns of Figure 26 if 
both ground and excited vibe-rotor states have similar RES-shape. However, only sub-class ii(180°) 
patterns (with ii < 0) exhibit spectral ordering {AiTiE){T2Ti){ET2A2){T2Ti) on the left hand side 
of Figure 26 that is maintained even as levels re-cluster into patterns {TiET2){TiET2){A2T2TiAi) 
of C3[iii] local symmetry across the separatrix break on the right-hand side of Figure 26 as analyzed 
below [8,37]. 0-crystal-field wavefunctions for either case tend to follow a Bohr-orbital progression 
s{Ai),p{Ti), d{E, T2), /(Ti, A2, T2),g{E, Ti, T2, Ai), ... In general, ordering is sensitive to RES-shape 
and tensor rank as discussed later. 



Table 11. Splittings of O D C4 given sub-class structure. 
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Figure 25. O i-class level clusters of C4 local symmetry (a) O4; (b) I4; (c) 24; (d) 84. 




For an isolated three-level (ATE") -cluster of local symmetry O4 or else 24 the splitting pattern requires 
only two parameters. This could be either the 180°(ii,iii) or the 90° (i?^;^,/?^) class pair in Table 11. The 
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120°-class, lacking 180° flips, has just one real parameter rj. Parameters zj, Rxy, and n each split (ATE) 
by 2:1 ratio but differ in sign. 

Local symmetry I4 and 84 each have two-level (TT) clusters that require just one splitting parameter, 
say ii, or else Rxy. Complex parameters and of the 90° i?^-class and the p^(180°)-class in Table 1 1 
may play minor roles in most C4 clusters but are necessary in order that the whole set be orthonormal 
and complete. 



Figure 26. Excerpts of SFq z/4P(88) superfine spectral cluster structure in 16/im region 
(Missing: K4=82...t 
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7.1.4. Orthogonality-Completeness of Local Symmetry Parameters 

Equation (122) expands P^n^ by Equation (83) in group operators (l,ri,r2, ...ie). It acts on |1) to 
give \iZn4) eigenkets in Equation (129). 



on / fa '± 

& = P^rJ 11) V = V ^ 5] D^nr^*i9>^)Sl, |1) 

6=1 



(129) 



6=1 
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An 0-symmetric H matrix is a sum of dual operators (1, Fi, f2, ...ie) with coefficients ga=^o, ri,r2, ...,i6. 
Local symmetry C4 or C3 reduces the sum to pg=10 sub-class terms Ca=ga+g[,+... each sharing a 
coefficient ga=g'a- 

°G PG 

H = Qaga = 9aCa (130) 

a=l a=l 



From these arise expansions like Table 11 of H eigenvalues e^^ in terms of its coefficients ga. Dual 
commutation gjgk=SkSj makes and H commute. Duality relation in Equation (94) leads to a 
L)^* -weighted sum of ga analogous to sum in Equation (129) of\ga). 

° G ° G 

e: = ei H It^) = (1| Pt)HP("J |1) ^ = (1| HP(<^J |1) ^ = {llY^gaga E Diti9t)gt |1> 

= {l\Y.9aT.^nJ*i9b)gbg-' |1) = = ^^^^(^a) °C„^7„ 

a=0 b=0 a=0 a=l 

(131) 

Each C4 sub-class of order has equal terms gaD^'^{ga) = g'^D^^n'' {g'a) =.. .expanding 
eigenvalue e^^. Rank-of-group = 10 is the number of eigenvalues and of expansion terms 
°CagaD^''{ga) in Equation (131) or Table 11. Each of ten eigenvalues e^^=(e^^6^2^ •••5^34) expand 
to ten C4-local tunneling parameters ga={^o, ^ih ••"> ki) and vice-versa. 



5. = (1| H = (1| Hg, |1) = E E E ^fkiOa) (1| HPJ, |1) 

a J k 

= E E {9a) (1 1 HPl 1 1) = E E (^?a) el 



One might count twelve real parameters in Table 1 1 since both pairs (ri,fi) and (Rz,Rz) are complex, 
unlike Ri = Ri, which are real. If i7 is a Hermitian array (H = it should only require ten, the rank 
of O, for its ten distinct real eigenvalues and the parameter pairs must be complex conjugates. 

With no conjugation symmetry, such as for a unitary O D C4-symmetric matrix, the R and r 
parameters may be complex and unrelated to R and f, and resulting extra real parameters are then 
needed. Symmetry parameter dimension matches eigensolution dimension for each local symmetry as 
shown in Figure 27. 

7.1.5. Resolving Hamiltonians with C3 Local Symmetry 

The previous two sections have detailed of symmetry-based level clustering and cluster splitting for 
C4. In Figure 26 these are the lower energy clusters of SFq for z/4 P(88). Given the previous two 
sections, it is possible to find the splittings of the Cs sub-group quickly. Starting with Equation (117) 
and Equation (118) one can build the irreducible representations necessary to create the Pnsns Ihe 
new sub-group. At this point, one can create a table analogous to Table 11. Such a table for Cs is shown 
in Table 12. The C3 clustering fits patterns of (^41, A2, T2, T2) and two of Ti, T2), each with a total 
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degeneracy of 8. As before in Figure 25, the splittings in C3 make different patterns depending on which 
tunneling parameters are active. This is demonstrated in Figure 28. 

Figure 27. OcL-local symmetry eigenmatrix parameters (a-e) L=C4,...,Ci (f-j) 
L^0,D^,...,D2. 
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Table 12. Splittings of O D C3 given sub-class structure. 
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Figure 28. O i-class and p-class level clusters of Cs local symmetry given different 
tunneling parameters. 




7.1.6. Octahedral Splitting for a Range of Local Symmetry C1CC2...CO 

As the order °L of local symmetry LcG decreases there are proportionally fewer types of local 
symmetry irrep d^{L) and hence fewer types of energy level cluster since each cluster is defined by its 
induced representation (i^(L)tG. There is a proportional increase in total number i^^^={£^yGrL of 
levels in each eigenvalue cluster. However, G-symmetry degeneracy limits the total number of distinct 
eigenvalues from all clusters to be global rank p{G) or less, no matter what local symmetry is in effect. 
Octahedral rank is p{O)=10=i^^-\-i'^'^-\-£^-\-£^^-\-£^'^ where £^ gives both the global degeneracy of each 
level type and the number of times it appears. 

The number of i7-matrix parameters equals the number of distinct eigenvalues as long as all 
eigen^'ectors are determined by global-local symmetry, that is, each entry is 0 or 1 in the GdL 
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correlation array. Diagonal eigenmatrix forms are shown in Figure 27a,b for C4CO and CsCO for 
which all bases states are distinctly labeled. Multiple correlation (> 2) occurs if L-symmetry is too 
small to determine some of the °G eigenbases. Then the H -matrix must have extra parameters that fix 
vectors through diagonalization. 

This happens for the C2(ii) cO symmetry whose correlation array in Equation (117) assigns the same 
C2 label to two bases of Ti and of T2. (Two C2 symmetries O2 and I2 cannot distinctly label three bases.) 
Figure 22 shows C2(ii) splits O into fourteen sub-classes: (1), (rif4), (r2f2), (rafa), (r4fi), (pxPz), 
(Py), (Ra,R^), (Ra,R^), (Ry'Ry)Xh)^ (12), (1315), (Uie)- The C2CO sub-classes form a non-commutative 
algebra and cannot be resolved so easily as C3CO or C4CO into commuting idempotent combinations 
like Equation (123). 

Spectral resolution of fourteen C2(ii)cO sub-classes requires more than rank number p(O)=10 of 
diagonal commuting O idempotents P^w To fully determine C2 basis, two off-diagonal pairs Pj^=P^^^ 
and Pj2=p^2t non-commuting nilpotent projectors are needed to finish C2-labeling of T-triplets. 
Adding these four gives fourteen projectors with their fourteen parameter coefficients shown in 
Figure 27c to fully define general C2(ii)cO i7-operators. (However, only twelve of the fourteen 
parameters are independent for Hermitian 77^,6=^5* «•) 

The other class of C2 symmetry has similar problems. Local C2(pz)cO symmetry requires projector 
pairs Pj^=P^^^ and P^l=P^^\ too, but then another nilpotent pair P^^=P^^ must be added to label 
repeated E bases in array Equation (117). This gives sixteen C2{pz) sub-classes to resolve and sixteen 
parameters sketched in Figure 27d. (Hermitian H=H^ matrices for C2(pz)cO have thirteen free 
parameters.) 

For the lowest local symmetry Ci=[l] (i.e., no local symmetry) sub-classes are completely split 
since every 0-operator is invariant to 1 as Ci provides no distinguishing labeling, and all twenty-four 
0-projectors (Sc(^'^)^=24) are active in its resolution. The 24-parameter i7-matrix resolution is sketched 
in Figure 27e. Each parameter for a=l, 24 is a combination of 24 products Dj'l{gp)gp (p=l, 24) 
of irrep and group element coefficient gp as given in Equation (90) or (131). (If i7 is Hermitian the 
number of free parameters reduces to T>a£^{i^+l)=n.) 

For O's highest local symmetry, namely O itself, there is no splitting of the T>a{£^)^=5 invariant 
idempotents P^ that resolve the five O classes. Then H has five independent parameters and five 
eigenvalues of degeneracy {£^y. This 5-parameter resolution is sketched in Figure 27f. Total level 
degeneracy for sub-matrix eigenvalues are listed below each one, and show less splitting than Abelian 
cases listed in Figure 27a-e. 

Any non-Abelian local symmetry such as L = D4 also fails to split P^ into a full number £^ of 
components P^^ if O irrep- (a) correlates with multi-dimensional L-irreps. By splitting out less than the 
full rank number p(O)=10 of idempotent projectors P^^, the resulting number of independent H matrix 
parameters reduces accordingly. The 8-parameter resolution for an i/-matrix with i^4C0 is sketched in 
Figure 27g and similarly for D^cO in Figure 27h. Two kinds of i^2C0 in Figure 27i,j share degeneracy 
sums with the Abelian cases. 

Each matrix display lists exact degeneracy due to global symmetry O but not the cluster 
quasi-dtgtntracy £^^^ due to local symmetry induced representation d^{L)'\'G. The latter is found 
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by summing global degeneracy £^ of all states |^;^) with the same local symmetry A as per Frobenius 
reciprocity in Equation (106). The result is integer ^^^^=(^^)°G/°L mentioned above. 

8. Spectral Resolution of full Oh Symmetry 

Including inversion I and reflection operations allows parity correlations between ewen-g (gerade) 
and odd-t^ (ungerade) states. Two classes of C2 subgroups lie in O and appear in separate C2 -correlations 
in Equation (1 17). In the following Oh correlations Equation (133), the two types of C2V subgroups have 
separate tables. The first subgroup C2^=[l, ay, ii, 0-2] is the one of the three local symmetries shown in 
Figure 12 while the second C|^=[l, p^, ay, ax] is just a subgroup of local symmetry C^y as would be 
Cl^=[l,p^,as,a4\. 
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The local symmetry C^^ydOh unambiguously defines all states in its correlation array while the other C2V 
symmetries fail to split the Eg and Eu sub-species. The former lead to complete eigenvalue formulae. 
The latter may not. 

8.1. Resolving Hamiltonians with C2V Local Symmetry 

As the order of the local sub-group symmetry goes down, the degeneracy and complexity of the 
rotational cluster must increase. Oh^C2v clusters are 12 fold degenerate and come in 4 cluster species. 
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Matrices describing this system are larger, but ODC2 will show many of the same effects. To actually 
resolve the doubled Ti or T2 triplets of requires distinguishing the u and g versions of each. The 

C2 clusters are 12 fold degenerate, but they are also easily displayed. 

As noted earlier, OdDsDC2 and ODD4DC2 local symmetries give identical cluster degeneracies and 
groupings, but with cluster splittings and structure dependent on the sub-group chain. Though it neglects 
inversion. Figure 27 indicates that there are several different types of ODC2 (and, thus Oh^C2v local 
sub-group symmetries). Examples given here involve the O D 04^02(14) sub-group chain. 

Compared with ODC4 and ODC3, the splittings of ODC2 are relatively simple to calculate since the 
terms in Equation (131) will be real. Creating splitting tables for C2 is done in the same way as for 
Tables 11 and 12. It is shown in Table 13. 



Table 13. Splittings of O D 02(14) given sub-class structure. 
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8.1.1. Local Sub-Group Tunneling Matrices and Their Inverse 

Table 13 can be further broken apart to demonstrate how one can create an automated process 
to evaluate the tunneling splittings for ODC2 local- symmetry structures. What will result is 
a transformation between cluster- splitting energy and tunneling parameters. The inverse of this 
transformation is also easily defined. 

Equation (131) produces Table 13, but even after combining splittings from each subclass, repetition 
exists. We show the two steps to convert Table 13 into the transformation matrix just described. First we 
assume that only levels may interact with themselves, e.g., that a I2 cluster may not interact with a 
O2 cluster. Second we recognize that only half of the subclasses are needed to fully define the possible 
splittings, the others simply repeat the same information. Table 13 shows this for the O2 cluster. Looking 
at the Ai level in the O2 cluster, one can see that the subclasses 1, r^, make a vector {1, 4, 4, 2, 1} 
while the i?^, in subclasses make a vector {4, 2, 4, 1, 1}. These vectors are reordered versions of each 
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other. Thus only one is needed. The A2 level in the I2 cluster shows the same similarity, but the i?^, in 
now contain a negative sign. 

By using only the minimum number of splitting parameters and including only a single cluster gives a 
condensed version of Table 13 that acts as a transformation that inputs symmetry-based tunneling values 
and outputs energy levels. Such a table is shown in Table 14. A simple inverse of the matrix in Table 14 
will produce the transformation giving tunneling parameters for a given set of cluster energy splittings, 
as shown in Table 15. 

Table 14. Matrix that converts tunneling strengths to cluster splitting energies. 
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Table 15. Matrix that converts cluster splitting energies to tunneling strengths. 
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There are multiple ways to use Tables 14 and 15. Among the most useful is to use the columns of 
Table 14 as a predictor of possible splitting patterns. Using the inverse matrix to find spectroscopic 
tunneling parameters from cluster splittings may also become a useful and automated process. 

An example demonstrates this process for a model (4, 6)-octahedral-Hecht spherical-top Hamiltonian 
Equation (134) with varying spectroscopic parameters. The terms and model rotational 
distortions written in an octahedral basis of fourth and sixth order respectively in J. The parameter 
9 is varied to explore the different relative contributions of and while keeping them normalized. 
Because and T^^l each have octahedral symmetry. Equation (134) represents all possible octahedral 
pure rotational Hamiltonians up to sixth order. 

H = BJ^ + cos(e)TW + sin(e)T[^l (134) 

As noted in Section 3 cluster structure location and the RES shape will change significantly as the 
Hamiltonian parameters change in Equation (134) as Figure 29 (a copy of Figure 6) shows by plotting 
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rotational energy levels of Equation (134) for changing 9 with corresponding RES at points along the 9 
axis. RES plots in the figure demonstrate how the phase-space changes as 9 varies. 

RES diagrams in Figure 29 along with the cluster degeneracy indicate where in the parameter- space 
C2 clusters exist. The lowest 02(C2)tO cluster in Figure 29 for 9 between 18° and 132° labels a 
kaleidoscope of 12 waves each with C2 local symmetry. Its superfine levels are magnified about 
100 times in the central inside plot of Figure 30 which has been adjusted to show level splittings 
but not whole cluster shifting. (The 6^-dependent cluster center-of-energy is subtracted.) The locally 
antisymmetric l2(C'2)tO clusters contain quite similar superfine structure but with A2 replacing Ai and 
Ti switched withr2. 

Figure 29. J=30 Energy levels and RES plots for T^^^^^vs.[4,6] mix-angle 9 with 
levels above 0=0° (extreme left), Tl^l levels at 0=90° (center), and -T^ levels at 6^=180° 
(extreme right). C4 local symmetry and 6-fold level clusters dominate at 0=17° while Cs 
type 8-fold level clusters dominate at 0=132°. In between these extremes are C2 type 12-fold 
level clusters particularly around 0=80° where a C3 — C4 level-cluster-crossing of the top 
14 levels occurs. 




At certain 0-points in Figure 30 levels of different symmetry cross and one of three distinctive 
splitting patterns emerge. These points occur periodically as indicated by vertical lines that are 
(starting form left side) solid, dotted, dashed, dotted, solid, dotted, dashed, solid, and so forth across the 
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plot. The three distinctive e^-energy level patterns for species a={Ai, E, Ti, T2, T2) are given by vectors 
Qa5/i=(0, 0, 0, 1, -1), erfot=(2, -1, 1, 1, -1) and esoUd^i"^^ -1, 1, 0, -1), respectively. These repetitious 
patterns seem to persist even outside of the marked-off sections to the very ends of the C2 cluster region 
at 0c^l8° and 6^c^l32° where they grow slightly but maintain their respective superfine ratio patterns and 
degeneracy. The matrix in Table 15 transforms each of the three e^-vectors in Figure 30 into a vector 
of 0-defined sub-class tunneling amplitudes g^. These are evaluated for clusters at several values of 
parameter 9 used in T'^'^l Hamiltonian Equation (134). Proportioned values of the tunneling amplitudes 
for the three distinctive cases are listed in the inset legend of Figure 30 based on Tables 14 and 15. 

Figure 30. The plot focuses on the lowest 02(^^2) tO cluster in the previous energy plot 
(Figure 29) of the Hamiltonian for J = 30. The inside plot has been magnified 

100 times. The inside diagram also centers the levels around their center-of-energy, showing 
only the splittings and ignoring the shifts of the cluster. Symmetry species are colored as 
before: Ai\ red, A2: orange, E2 : green, Ti\ dark blue, and T2: light blue. The vertical lines 
on inside plot draw attention to specific clustering patterns described in the text. l2(C'2)tO 
clusters have similar superfine structure but with A2 replacing Ai and Ti switched with T2. 
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Dotted-line and solid-line curve patterns appear alternately flipped in sign. Dotted-line patterns have 
a crossing (Ti,T2) pair while the solid-line patterns have a crossing (T2,E) pair. 

Solid-line patterns appear to be centered on quasi-hyperbolic avoided-level-crossing episodes 
involving the pair of repeated T2 tensor species of O. The ordering (^41, Ti, T2, T2) of solid-line 
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superfine level patterns reflects Bohr-like orbital ordering {s,p,d, f , ..) of orbital momentum and occurs 
only when there is just one non-zero sub-class of tunneling parameter, namely that of sub-class (r34 or 
equivalent Rxy) that affects tunneling between nearest-neighbor C2 valleys. 

Dashed-line pattern level curve slopes appear to alternate (+) and (-) signs and exhibit maximum 
separation of repeated r2-species surrounding a degenerate (yli, Ti, £") -sextet crossing midway in 
between. Such triple-point crossings are quite remarkable. They appear repeatedly in Figure 30 and 
persist even at low- J as seen for J=4 in Figure 8c. Higher 12(^^2) tO clusters show similar triple points 
made of {A2, T2, £") -sextets. 

Such crossings are quite ironic if we recall that it was {AiTiE), (A2T2E), and (T1T2) clusters 
noted by Lea, Leask, and Wolf [20] and later Dorney and Watson [21] that led to a theory 
involving induced representations K4(C4)tO including 04{C4)'tO=Ai®Ti®E, 24{C4)'tO=A2®T2®E, 
and ±l4(C4)tO=Ti0r2. (Recall C4 columns of Equation (117) and reciprocity Equation (106)). This 
theory uses an inter-C4-axial tunneling model [22,23] with a single ad. hoc. tunneling parameter that 
predicts a 2:1 -splitting ratio for (ATE) clusters. C4-axial tunneling cluster splitting dies exponentially 
as body momentum-X approaches J (Recall Figure 26) and thus C4{ATE) levels never actually cross. 

However, the C2{ATE) levels in Figure 30 clearly do so and with quite the opposite l:2-spltting ratio. 
It is ironic that the more elegant ortho-complete multi-path tunneling models, while useful in exposing 
these crossings, seem at a loss to explain them, particularly given that they were first noted by Lea, 
Leask, and Wolf so very long ago! 

It would be easy to write off such (ATE) triple-crossings and particularly the (T1T2) or (ET) 
double-crossings as "accidental" degeneracy. Indeed, all but the latter occur for special values of a 
complete set of sub-class parameters. However, Figure 30 clearly shows that each type of crossing 
belong to a periodic structure that is unlikely to be just an accident. 

Clearly there is still much to learn about multi-path tunneling models in general and the octahedral 
ones in particular. Here we can only offer a potentially elegant way to treat these kinds of 
high- symmetry cases. 

9. Examples of Rovibronic Energy Eigenvalue Surfaces (REES) and J-Clusters 

Semiclassical treatment of rovibronic or rovibrational states provides some insight into the transition 
between lab-coupled and body-coupled vibronic momentum that are related in Equation (8) through 
Equation (10b) of Sections 1 and 2. The first semiclassical analysis of fundamental coupling in high- J 
octahedral molecules was done for iy2E [38] and [39] bands in 1978 and for overtone 1^2 + ^3 
"hot-bands" in 1979 [40]. 

These methods are similar in philosophy to those described in Section 2 that approximate tensor 
eigenvalues with Legendre formulas and thereby construct rotational energy based on a semiclassical 
J-vector. However, the more general approach differs in that it builds an A^-by-A^ matrix of such 
formulas that takes account of quantum rovibronic coupling between vibronic (or vibrational) states, 
that is, a 2-by-2 matrix for the U2E system, a 3-by-3 matrix for the ^^3X1 system, and a 5-by-5 matrix for 
the U2 + ^3 system. 
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The resulting eigenvalues provide points on nested Rovibronic Energy Eigenvalue Surfaces 
(REES) for each direction of the semiclassical J-vector. Visualization of P, Q, and R state mixing in 
bands by 3-sheet REES was done using the high-resolution 3D-graphics at Los Alamos in 1987 and 
reported in 1988 [25]. Interesting features of the Z/3T1 REES include conical intersections that occur for 
zero scalar Coriolis coupling. These are analogous to well known conical intersections of Jahn-Teller 
PES that lend insight into BOA breakdown of single adiabatic surfaces. The following contains two 
examples of REES models. The first is a simplified internal rotation model involving a 2-sheet REES, 
and the second is an excerpt of a recent study of the pz/2tj^ dyad of CF4 that involves a 9-sheet REES. 

9.7. Rotor-With-Gyro Model of Internal Rotation 

A first application by Ortigosa and Hougen [17] of REES to visualize molecules with internal rotation 
is related to a simple rotor- with-gyro model [25,41] based on the three lowest rank tensors possible, 
namely the scalar (rank-0), the vector (rank-1), and the tensor (rank-2). The prolate symmetric top RES 
in Figure 1 is an example of a scalar-tensor combination. A vector RES lacks J-inversion symmetry, 
that is, time reversal symmetry, so it is forbidden for normal molecules that have no intrinsic dynamic 
chirality such as embedded spin S. We consider how to include an S in a way that preserves overall 
T symmetry. 

Total momentum J=R+S is the sum of rotor momentum R and gyro spin S. J is conserved in lab 
frame but R and S are not. If gryo is body-frame-fixed by frictionless bearing then rotor gyro-coupling 
does no work and is an ignorably constant Hrs- S and |J| are conserved in body frame but J and R 
are not. 



Replacing bare-rotor momentum R=J-S gives the following with a new constant spin energy H' rs. 



The simplest classical theory of the rotor- R-gyro-S momentum dynamics involves superimposed RES 
plots, one for +S and one for -S in Figure 31; A composite RES with T symmetry. If J and +S align 
(anti-align) then |R|=|J-S|, rotor energy Equation (135), and rotor-gyro relative velocity are 
minimized (maximized) (Thus, gyro-compass alignment with Earth rotation is seen to be relativistic 
quantum effect!). 

A quantum theory of multiple RES involves mixing extreme cases | J ± S|. An elementary quantum 
gyro-spin is a two-state spin- 1/2 with a 2-by-2 Hamiltonian matrix found by inserting quantum spin 
S=(7/2 matrices into Equation (136) to give Equation (137). Gyro-rotor dynamics involves REES 
obtained from eigensolutions of the following 2-by-2 matrix for each body-based J-vector Euler 
orientation 7). 



Hrj^S {hod- fixed) — 



+ BY{? + CR^ + Hs + Hrs 



(135) 




(136) 
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Hr^S (quantized) = + 5jJ + CJ^ - Ai^G^ - BiyGy - CJ^G^ + COUSt. 

(RErotor — JC COS /3 —AJ COS 7 sin /3 — iBJ sin 7 sin /3 

—AJ cos 7 sin ;5 + i5 J sin 7 sin ;5 RErotor + JC cos /5 

where: RErotor = J^(^cos^7sin^/3 + 5sin^7sin^/3 + Ccos^f3) 

Eigensolutions of matrix form Equation (137) transform classical RES Figure 31 into quantum REES 
Figure 32 that has conical intersections or avoided crossing points replacing lines of classical surface 
intersections in the former Figure 31. Also, individual sheets of REES have J-inversion symmetry (or T 
symmetry) that individual RES lack. Where the RES of Figure 3 1 are well separated their shape is not 
so different from that of REES in Figure 32. Differences show up near the intersection lines where the 
two RES approach resonance. In this resonance region the REES is deformed extremely from rank-1 or 
rank-2 tensor shape of the separate RES, and there arises greater mixing of the extreme |±S| base-states. 

Figure 31. Views of classical rotor-gyro RES for spin +S (yellow) and — S (gray). 





Figure 32. Same views of quantum REES for rotor with gyro spin operator S=g/2. 
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9.2. REES ofCF^ in v^/2v^ Dyad 

The first practical REES application includes 9-sheet displays of the pz/2v^ dyad of CF^ recently 
shown by Boudon etal. [2]. This large scale numerical analysis may be summarized by a revealing plot 
of dyad eigenlevels as a function of J = 0 to 70 in Figure 33. This includes colored lines representing the 
REES values for J located on C4 axes (shaded red), C3 axes (shaded blue), or C2 axes (shaded green). 

Each of the symmetry axes may take turns as central loci for clusters of their type of local symmetry 
^2,^3, or C4, or else, they may sit on a REES separatrix or saddle point between two or more different 
types of clusters. A third option involves Ci clusters that have no rotation axis point but are likely 
to belong to vertical xyz-plant reflection symmetry Cy = [1, Pz] or diagonal-plane reflection symmetry 
Cd= [1, ^3]. These label clusters of 24 levels associated with 24 equivalent REES hills or valleys. 

Figure 33. (After Boudon et.al. [2].) (J<70) rotational levels of 1^3/2^4. 




J 

A final option involves true-Ci clusters with no local symmetry whatsoever and 48 REES hills or 
valleys. So far this extreme type has not been identified, but one may speculate that it may actually 
become most common at extremely high J. 

A common ordering noted before on the left hand side of Figure 29 (pure T'^l) and in Figure 26 
(16/x region of SFq) is (C3-valley^C2-saddle^C4-hill). It is present in the lowest REES band of 
Figure 33. An inverted version of the common ordering appears clearly in the 2^^ band whose REES is 
cubic in Figure 34. 

A cutaway view at J = 57 of the first five REES sheets shows glimpses of the first two REES deep 
inside of Figure 34. The second sheet has cubic topography similar to the inverted T'^l RES on the right 
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hand side of Figure 29 (pure (— )TW). However, the first and lowest REES for J = 57 is practically 
spherical with all2J+l=115 levels and clusters crushed in Figure 33 into near degeneracy! 

After the first two REES sheets the cluster topography become more complicated with multiple 
conical intersections and avoided crossing points. 

On the 5*^ sheet of the ( J = 57)^3/2^4 REES are found examples of Ci -local symmetry valleys as 
shown in Figure 34. (The upper four sheets are made invisible.) Each Ci loop occupies an area that 
is comparable to the minimum uncertainty (J = K)-cont shown on vertical C4 axis of the figure and a 
nascent 24-level cluster of type 12(^^2 )tO should be present in the level spectrum. 

The symmetry details in this rovibrational spectra and the potential richness of quantum dynamics it 
represents should be quite evident from the few examples glimpsed here. We seem to be just scratching 
the surface of quantum systems of a great but potentially comprehensible complexity. 

Figure 34. (After Boudon et.al. [2].) A rare (J=57)h{C2)"tO structure on fifth REES. 




10. Summary and Conclusions 

Semiclassical methods for visualizing and analyzing rovibrational dynamics of symmetric polyatomic 
molecules have been reviewed. This includes improved understanding of RES and REES phase spaces 
and development of more powerful symmetry methods to calculate tunneling dynamics of symmetric 
molecules that are highly resonant. A group-table-matrix analysis of intrinsic vs. extrinsic symmetry 
duality (The "Mock-Mach-Principles" Equation (95) and Equation (94) of wave relativity.) leads to 
generalizing character relations between group classes and irreducible representation into sub-character 
relations between sub-classes and induced representations Equation (131) and (132). These provide 
ortho-complete parameter relations (Tables 11-15) for complex tunneling path lattices that determine 
molecular fine, superfine, and hyperfine spectra. The methods may be extensible to fluxional atomic and 
molecular systems. 
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Appendix 

A. Classical Ds Modes: Local C2 and C3 Symmetry Examples 

Local symmetry theory applies to classical vibrational modes as well as to quantum tunneling. 
Examples of classical modes given below help clarify global-^5-local symmetry and geometry of 
group projection. For example, modes defined by local 02(^3) in Figure 20 are to be compared with 
quantum waves in Figure 19. Each mode ket Ij^.) has the same coefficients D^^{^) for projections in 
Equation (78) as the waves do, but the mode shapes clearly display a vector geometry. 

In particular, global x-vector modes and (left Ei column in figure) "point" along global 
x-direction while vector modes \y^) and \y^) (right Ei column) "point" along global y-direction. Each 
global pair = x,y) is projected to be an ia-symmetric-antisymmetric pair like lab unit 

vectors [|x),|^)] (Recall Equation (78)). 

\x) = P°3 |l) ^ = (|l) + lig)) 1^) = pi3 |i) ^ = _ |i3)) (138) 

This exposes easy derivations of £'-irrep D^^{^)={j\^\k) in Equation (89). Irreps in Equation (87) 
such as D^^{y) for 120°-rotation r simply contain direction cosines (j'|r|/c)=ej«e^./e of rotated vectors 
[r|x),r|^)] relative to original (Note transpose of equation array to matrix array). 

{x\v\x) {x\v\y) \^(-\ ^ 
{y\v\x) {y\v\y) ) ^ +^ -\ j 

(139) 

This also fixes local transformations. Local x-vector modes and |^^) (lower Ei row in 
figure) "point" along local x-axes that are local radial lines while local ^-vector modes and \y^) 
(upper El row) "point" along local y-axes that are local angular lines. If global symmetry meets local 
anti-symmetry as in (or vice-versa in \y^)), a zero appears on the i3-axis in Figure 20. Singlet 
modes |^^) and \y^) avoid such conflicts by being all one or the other. 

For group-defined cases like Figure 20, symmetry arguments alone determine normal modes that 
usually require diagonalizing a K-matrix (below) just as tunneling states (Figure 19) usually require 
diagonalizing an i/-matrix. 

A. 7. Comparing K -Matrix and H -Matrix Formulation 

Classical modes are eigenvectors of force-field matrix K or operator K that is a linear function of 
spring constants (/cq, etc, in Figure 35) for a harmonic approximate potential l^(x) that is a quadratic 
K-form of coordinates Xa based on six i^a-labeled axes or \a) shown in Figure 20. Each K component 
Kah={ci\i^\b) is a sum over spring fc-constants that connect axis-x^ to axis-x^ multiplied by factor 
(k^ ©x"^) (k^ ©x^) for projecting spring /c's end vectors k^ and k^ onto x"^ and x^ at respective connections. 
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(A straight- line spring has equal k<j=kft. Curvilinear springs must only have k-ends with equal sense 
(—)■—)■) or (■^■^) of spring direction. Either direction gives the same Kab)- 

V{x) = '^^^ {x\K\x) where: \x) = x"" \a) , (a, 6) = (1, r"^, r^, «i, i2, ^3) 



(fc) 



^^KabX^x^ where: Kab 



(fc) 

-EMk.»x«)(k,.x'') tf:aj^b 

(k) 



if : a = b 



0,6 



(140) 



This sum of harmonic Hooke (/ca;^/2) -potentials has diagonal Kaa terms followed by off-diagonal terms 

{Kab^ Kba). 



(k) 



2 

(fc) a,b 



k 



(k) 



(k) a=^b 



(141) 



Figure 35. X3 spring models with local symmetry: (a) 03302(13); (b) Mixed. 
(a) Local DpC2(ij} model (b) Mixed local symmetry model 

cosl5°=V[(l+cos30°)/2] 
=(I/2)V(2+V3)=sin75° 





sinl5°=V[(l-cos30°)/2] 
> -(l/2)V(2W3)=cos75° 



The classical equation of coupled harmonic motion is a Newtonian F= M a relation of a 
n-dimensional force vector F, acceleration vector a, and mass operator M. The latter is a 
unit-matrix-multiple M-1 for the i^a-symmetric case treated here. The driving force F is a (-)derivative 
of potential Equation (140) that becomes a K-matrix expression. 



dV 



(142) 



It is instructive to compare this classical equation of motion to that of Schrodinger's equation for 
quantum motion. 

i.b 



(143) 



Squaring quantum time generator ihdt=ii yields equations having classical form Equation (142) with 
K=H^ and M=h'^. 



-h'^d^i)'' = E Kabi'^ where: K = 



(144) 
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The (H//i)-eigenvalues are quantum angular frequencies em/h=^m' The (K/M)-eigenvalues are 
classical squared angular frequencies km/M=uj'^. The former is Planck's oscillator frequency relation 
e=lijjj. The latter is Hooke's relation k/M=(jj^. Apart from normalization, eigenvectors of quantum H are 
identical to those of classical K and either eigenvalue set corresponds to the respective energy spectrum. 

A.2. Comparing K-Matrix and H-Matrix Eigensolutions for Local DsDC2{is) 

The preceding relates eigensolutions Equations (92) and (93) of quantum Hamiltonian i7-matrix in 
(71) with those of a classical K-matrix. In particular, eigenvectors of H found using L^-matrices in 
Equation (89) or (139) also serve as mode-eigenkets in Figure 20 that diagonalize a 1^33(72(^3) -locally- 
symmetric X-matrix. With this symmetry, K cannot couple radial (local-x) and angular (local-y) modes 
and is left with just four independent real group-based parameters ga=ro, ri, i^, and is allowed for 
^3 ^^^2 (^3) -symmetric H in Equation (93). These relate to four spring fc/^-constants in Figure 35a. 

Only P^-row parameters 5'5=(l|K|g5)=Xi5 of the force matrix Kab are needed for the spring model 
in Figure 35a. That model includes /c^(angular) and /c^(radial) constants for internal connections between 
masses. The /c3(angular) and /co(radial) constants represent external connections between each mass and 
an outside lab frame. 

Generic group parameters gb=Hib^ labeled [ro, ri, r2, ii, ^2, ^3] for the i7-matrix in Equation (71), are 
now applied to gb=Kib. The gb are to be related to spring-constants kj using coordinate- spring projection 
cosine factors (ki • x^)(k5 • x^) in Equations (140) and (141). The usual harmonic limit assumes small 
vibrational amplitudes (x5<Cl) for which direction of spring end vectors or k^ do not vary to P^-order, 
and so, for lab-fixed the Kab are constants. 





ll> 




|r2) 


|ii> 


|i2) 


lis) 




ki/2 


ki/2 


ki/2 


ki/2 


ki/2 


ki/2 


(l|Kb) = 


-\-kj^ 


kf 1 2 


kj- 1 2 


+kr/2 


+kr/2 


krf 


+h 


+0 


+0 


+0 


+0 


-ks 




+ko/2 


+0 


+0 


+0 


+0 


+ko/2 



(145) 



One may visualize each —Kib as the acceleration of xi due to setting a (tiny) unit Xb in Equation (142). 
Diagonal -Ku must be negative or else xi blows up. Higher order anharmonic terms are needed 
to describe effects of rotating k^ or x^ and such models are likely to suffer from classical stochastic 
(chaotic) motion. 

Substitution of generic ga from Equation (145) into reduced £^33(72(^3) -symmetric i7-matrix in 
Equation (92) or (93) gives K-matrix eigenvalues K^^ due to each spring ki, kr, ks, or ko in Figure 35a 
separately or together. Modes in Figure 20 remain ei^^enmodes for all values of four spring constants ki, 
kr, ks, and ko since none can mix local x-and-y-symmetry. 

ro + ri + + + ^2 + ^3 = 
ro + ri + rl - ii - i2 - is = 





ri - - - ^2 + 2^3 - ri + - + ^2) 

■rl + ri - ii + ^2) 2ro - ri - + ii + ^2 - 2^3 
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Any set of four X-matrix eigenvalues /c^S /cf , and ky is arithmetically possible by adjusting the 
four spring constants. However, their arrangement in Figure 20 (this was drawn to match tunneling 
states in Figure 19) is impossible without negative fc- values that would give classical instability. As 
shown below, free ring molecules often have -stretching modes among the highest frequencies. In 
contrast, tunneling amplitudes are often negative so their Ai states lie low. As a rule, fewer quantum 
nodes imply lower energy. 

A. 3. K -Matrix Eigensolutions for Broken Local Symmetry 



In some ways the direct- fci -connection spring model of Figure 35b is quite the opposite of the 
D^'DC2{h) model just treated since it involves maximal (50-50) mixing of x and y local symmetry. 
Below are recalculated generic 5'5=(l|K|g5) in terms of direct spring-constants ki using (141) with 
projection cosines listed in Figure 35b. 





|1) 


|ri) 


|r2) 




|i2> 


lis) 


(1|K|^6) = 


A:i(cos2 75° 
+ cos2l5°) 
= A;i 


ki cos 75° 
• cos 15° 

_ fci 

4 


ki cos 15° 

• cos 75° 

_ fci 

4 


ki cos 15° 
■ cos 15° 
_ A:i(2- V3) 
4 


ki cos 75° 
■ cos 75° 
_ ki{2 + V^) 
4 


fci(cos2 75° 

- cos^ 15°) 
_ fci 

2 



(147) 

Again, a substitution of generic Qa from Equation (147) into reduced i7-matrix Equation (93) gives a 
reduced X-matrix like Equation (146), but now the E'-symmetry submatrix is not diagonal. 



1 

XX 



yy 

Kx ^ 



3A;i 
0 



xy 
E 



yy 



3fei 3fei 

4 4 

3fei 3fei 

4 4 



(148) 



Eigenvectors of the £'-submatrix are symmetric (+) and antisymmetic (— ) mixtures of x and y local 
symmetry states. 



K 



K 






(149) 



{x,y). 



Figure 36 shows (50-50 ±)-mixing due to ki. It distinguishes genuine vector modes (|f (_)) or |^(_))) 
and the scalar breathing mode (l^^^^.)) from non-genuine (low or zero-frequency) vector modes of pure x 
or ^/-translation (|f (+)) or |^(+))) and rigid rotation (pseudo-scalar |^^^)). The i3-local symmetry is wiped 
out by direct connection- fci. 

In order to reestablish approximate 1^33(72(^3) -local-symmetry there needs to be a (72(^3)- "/oca/e" 
provided by lab-grounded potential springs such as those with constants k^ and fco in Figure 35a. Adding 
these in the form of Equation (146) to Equation (148) causes a transition between the two extremes. If 
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the difference {ks + 2kr — ko) between eigenvalues Equation (146) begins to dominate the off -diagonal 
component (3/ci/4) of Equation (148), then mixed £'-modes of Figure 36 begin to recover 0^^02(13) 
locality seen in Figure 20 . 

Meanwhile the constant ks that determines eigenvalue fc^^ joes not affect locality for either of the 
singlet Ai or A2 modes. Singlet eigenvectors are non-negotiable as long as master symmetry holds. 

Figure 36. Mixed-local symmetry modes of direct- fci -coupled X3 model in Figure35b. 




A,4, K -Matrix Eigensolutions for D^dCs Symmetry 

Another choice for local symmetry is the C3 subgroup of Equation (79) corresponding to a strong 
chiral perturbation by internal rotation, spin, or 5-field. The £'-submatrix in Equation (146) with zero 
generic reflection parameters (^1=^2=^3=0) may take a purely chiral C3 form if the generic rotation 
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parameters ri and r2=r\ are purely imaginary corresponding to velocity dependent force ( ri=ir and 
r2=-ir. Here K is assumed Hermitian self-conjugate as was H). 



1 

XX 



yy 



1 XX 


xy 






^yy ) 





+ir\/3 



(150) 



Ti = ir = 
H = ^2 = ^3 = 0 



Cs ^'-eigenvectors have local x ± iy={±l)3 combinations that exhibit purely circular right R=(+l)^ 
and left L=(— 1)3 polarization orbits of C3 symmetry shown in Figure 37 (Recall C3 splitting in 
Equation (79)). 

Figure 37. L^a DCs-local symmetry modes of X3 molecule. 

'J 



Scalar mode 

A, 



Co coupling 




Genuine vibration mode^ 



Vector Right polarized 

^1 
R(-) 



Low-frequency modes 





Right circular-translation 



Psuedo- rotation 
A, 



Left circular-translation 
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K 



K 





(151) 



9 = (^,?/). 



Pure Cs symmetry is a normal subgroup and restricts fc^J, and fc^^ become degenerate. Both the 



scalar 



03,03/ 



and pseudoscalar 103^03) ^^^^^ t)Oth labeled equally by (0)3 symmetry. Local symmetry 
effectively goes global in the pure Ca-case where all internal coupling is zero. 

Any internal or external parameters may split the A1-A2 degeneracy and mix the C3 states to 
form elliptical polarization orbits. This is most efficiently calculated using U(2) analysis similar to 
Equation (137). 
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